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SECTION  1.  INTRODUCTION 


The  accuracy  of  a  free-flight  rocket  depends  significantly  on  the 
transverse  angular  rate  of  the  rocket  at  the  time  it  is  no  longer 
physically  constrained  by  the  launcher  guidance  mechanism,  i.e.,  at 
end  of  guidance  (EOG) .  If  the  causes  of  such  angular  rates  cannot  be 
accurately  predicted  and  proper  allowances  made  for  them,  an  accurate 
free-flight  rocket  system  cannot  be  constructed,  except  perhaps  by 
trial  and  error.  Historically,  launcher  motion  and  guidance  mechanism 
imperfections  have  been  considered  principal  causes  of  transverse  angular 
rates  of  free-flight  rockets  at  EOG.  Dynamic  imbalance  of  spinning 
rockets  and  thrust  misalignment  have  also  been  recognized  as  resulting 
in  what  are  effectively  secular  transverse  angular  rates  at  EOG. 

Due  to  difficulties  in  interpreting  optical  lever  data  in  the  light 
of  only  the  recognized  error  sources  mentioned  above,  and  to  the  fact 
that  the  data  in  question  evidenced  considerable  transverse  bending  of 
rockets,  both  during  and  subsequent  to  guidance,1  it  was  suspected  that 
there  was  a  strong  correlation  between  transverse  vibration  of  a 
free-flight  rocket  during  guidance  and  its  transverse  angular  rate 
subsequent  to  EOG.  Through  a  previous  research  effort,2  it  was  shown, 
on  the  basis  of  a  simple  two-body  model  of  a  spinning  flexible  rocket, 
that  if  such  a  rocket  is  vibrating  transversely  at  EOG  then,  in  general, 
it  will  possess  a  non-zero  transverse  angular  rate  subsequent  to  EOG. 
Furthermore,  two-body  model  results  indicate  that  the  transverse  angular 
rate  attributable  solely  to  transverse  vibration  may  be  on  the  order  of 
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100  mrad/sec.  To  put  this  in  perspective,  it  is  desirable  that  the  initial 
transverse  angular  rate  of  a  free-flight  rocket  (which  cannot  be  treated 
as  a  bias)  be  no  greater  than,  say  10  mrad/sec.  It  follows  that  the 
problem  of  transverse  vibrations  must  be  dealt  with  if  such  a  system  is 
to  be  achieved. 

Recent  research  has  provided  a  more  detailed  model  (assumed-modes 
model)  of  a  spinning  flexible  rocket  which  has  been  used3  to  validate  the 
two-body  model  and  also  to  obtain  certain  results  presented  herein.  In 
addition  to  the  development  of  the  aforesaid  model,  the  following  objec¬ 
tives  were  established  at  the  beginning  of  this  effort:  (1)  to  gain  a 
better  understanding  of  the  mechanism  whereby  transverse  vibrations  are 
transformed  into  transverse  angular  rates,  (2)  to  determine  probable 
causes  of  transverse  vibration  and  (3)  to  recommend  ways  of  reducing  trans¬ 
verse  vibration  during  guidance  and/or  its  adverse  effects.  The  following 
sections  of  this  report  describe  the  approaches  taken  in  pursuit  of  these 
objectives  and  the  principal  results  of  this  study. 

In  Section  2,  the  effects  of  transverse  vibration  which  is  present 
at  EOG  are  considered  without  regard  for  the  cause,  or  causes,  of  such 
vibration.  The  two  models  mentioned  above  are  described,  and  physical 
explanations  for  the  post-launch  transverse  angular  rate  generally  associated 
with  transverse  vibration  are  given.  Some  results  obtained  using  the 
mathematical  models  based  on  the  two-body  and  assumed-modes  models  are 
presented.  Methods  for  analyzing  optical  lever  data  not  considering 
vibration,  and  also  considering  it,  are  described  in  Section  3.  Probable 
causes  of  transverse  vibration  of  a  free-flight  rocket  during  guidance  are 
discussed  in  Section  4.  Results  are  presented  for  two  probable  causes, 
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support  imperfections  and  friction.  Conclusions  and  recommendations 
follow  in  Section  5.  Most  of  the  mathematical  details  are  relegated  to 
appendices . 

The  term  "mallaunch"  is  avoided  in  this  report,  because  the  meanings 
attributed  to  it  by  various  individuals  are  not  the  same.  Instead,  we 
shall  use  the  expression,  "transverse  angular  rate  at  end  of  guidance" 
in  which  "transverse  angular  rate"  refers  to  the  magnitude  of  the  angular 
velocity  of  the  instantaneous  centroidal  principal  axes  of  the  vibrating 
rocket  and  "end  of  guidance"  is  an  event  which  occurs  at  the  time  of  last 
intended  physical  contact  between  the  rocket  and  its  launcher. 


SECTION  2.  EFFECTS  OF  TRANSVERSE  VIBRATION 


2.1  Description  of  Models 

As  stated  previously,  two  models  of  a  spinning  flexible  rocket  have 
been  developed  and  used  to  study  the  effects  of  transverse  vibration  of 
such  a  rocket  during  the  guidance  phase  of  its  flight  on  its  motion 
subsequent  to  EOG.  One,  the  two-body  model,  is  fairly  simple  -  as  simple 
as  possible,  it  is  felt,  if  the  effects  of  transverse  vibration  are  to  be 
modeled.  The  other,  which  we  call  the  continuous  model  and  also  the  assumed- 
modes  model,  is  more  detailed. 

The  Two-Body  Model 

The  two-body  model  consists,  as  its  name  implies,  of  two  "cylindrical" 
rigid  bodies  (see  Fig.  1)  which  are  connected  by  torsional  springs,  viscous 
dampers  and  a  "hinge".  Then  allow  "transverse"  angular  rotation  of  the 
bodies  relative  to  each  other.  Here,  "transverse"  is  used  to  indicate 
motion  orthogonal  to  the  longest  axis  of  either  body.  In  most  of  the  work 
which  has  been  done,  both  bodies  have  been  assumed  to  be  axisymmetric  and 
homogeneous.  These  assumptions  could  have  been  avoided,  but  to  do  so 
would  have  resulted  in  a  more  complex  mathematical  model,  and  they  are 
very  proper  for  the  type  of  rockets  considered. 

As  shown  in  Fig.  1,  the  aft  body  is  denoted  by  Body  1  and  the  other  by 
Body  2.  Body  j's  center  of  mass  is  C j .  The  hinge  point  is  R  and  the 
points  P  and  Q  are  points  on  the  geometric  centerlines  of  Bodies  1  and  2, 
respectively.  The  points  P  and  Q  are  constrained  to  move  rectilinearly 
while  the  rocket  is  on  the  launcher,  if  the  mechanism  which  supports  the 
rocket  is  rigid.  If  the  supporting  mechanism  is  flexible,  as  is  assumed 
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to  obtain  some  of  the  results  given  in  Section  4,  then  the  points  P  and 
Q  are  constrained  in  another  manner  (see  Appendix  D) . 


Fig.  1.  Two-body  physical  model,  (a)  On  launcher,  (b)  Off  launcher. 

The  bodies  are  truly  rigid  in  that  their  masses  are  constant.  Since 
the  time  period  considered  is  that  from  ignition  to  shortly  after  EOG, 
this  assumption  is  a  good  one  if  the  rocket's  mass  does  not  change  substan¬ 
tially  during  this  short  time  period.  The  constant  mass  assumption  is  also  necessary  if 
the  closed-form  solution  to  the  equations  of  rotational  motion  given  in 
Appendix  A  is  to  be  obtained.  Although  the  mass  is  assumed  constant,  a 
constant  magnitude  thrust  force  (F^)  is  assumed  present.  The  thrust 
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acts  through  point  P  and  is  not  necessarily  aligned  with  the  longitudinal 
axis  of  Body  1.  That  is,  mechanical  thrust  misalignment  may  be  modeled. 
Another  basic  assumption  which  allows  the  closed-form  solution  to  be 
obtained  is  that  aerodynamic  reactions  (force  and  moment)  are  negligible 
during  the  time  of  interest.  This  assumption  is  a  good  one  if  the 
magnitude  of  the  velocity  of  the  rocket's  center  of  mass  is  "sufficiently" 
small  in  magnitude,  say  less  than  100  m/ sec . 

The  spin  rate  of  the  rocket  is  also  assumed  constant  to  obtain  the 
closed-form  solution,  Eut  it  is  modeled  as  variable  when  the  effects  of 
friction  (Section  4  and  Appendix  D)  are  considered. 

The  Assumed-Modes  Model 

The  name  "assumed-modes  model"  was  coined  because  certain  mode  shapes 
were  "assumed"  in  mathematically  modeling  (see  Appendix  B)  the  transverse 
vibrations  of  the  rocket,  which  is  modeled  physically  as  a  continuous 
flexible  body.  Two  sketches  of  the  continuous  physical  model  are  shown  in 
Fig.  2.  Figure  2a  shows  how  during  the  guidance  phase  two  points,  P  and 
Q,  on  the  deformed  centerline  of  the  rocket  model  are  constrained  to  move 
rectilinearly.  After  EOG  (see  Fig.  2b)  the  points  P  and  Q  may  move 
arbitrarily.  The  rocket's  mass  is  assumed  to  be  constant  and  a  thrust 
force  (possibly  time  varying)  acts  on  the  rocket.  The  axial  component  of 
this  force  is  assumed  to  act  at  the  head  end  of  the  rocket  motor  section; 
however,  the  effects  of  transverse  vibration  in  turning  the  flow  of  the 
rocket  exhaust  gases  as  they  exit  and  mechanical  misalignment  of  the 
rocket  nozzle  are  modeled  by  transverse  thrust  components  which  act  at 
point  P.  An  additional  characteristic  of  the  continuous,  or  assumed-modes, 
model  is  that  the  undeformed  shape  of  the  rocket  can  be  "nonstraight," 
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thereby  allowing  for  imperfections  in  motor  section  to  warhead  section 
connections  or  "non-straightness"  of  the  rocket  for  other  reasons. 


Fig.  2.  Continuous  physical  model,  (a)  On  launcher,  (b)  Off  launcher. 

The  transverse  vibration  of  the  continuous  rocket  model  during  guidance 
is  modeled  mathematically  using  pinned-pinned-free  normal  mode  shapes  of  a 
uniform  beam. 2  At  EOG,  a  transformation  is  made  so  that  free-free  bending 
mode  shapes  of  a  uniform  beam  can  be  used  for  the  free-flight  portion  of 
the  rocket’s  flight.2 
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In  the  assumed-modes  model,  the  spin  torque  may  be  time  varying. 

Thus,  a  spin  producing  device  such  as  an  eroding  spin  turbine4  can  be 
modeled . 

In  Section  4,  the  effects  of  supporting  the  rocket  at  three  points 
during  guidance  are  studied.  Although  this  could  be  done  by  using  pinned- 
pinned-pinned-free  normal  mode  shapes  for  a  uniform  beam,  the  approach 
taken  was  to  constrain  a  third  point  (between  P  and  Q)  to  follow  a  pre¬ 
determined  path  as  the  rocket  moves  longitudinally. 

2.2  Causes  of  Transverse  Angular  Rates 

There  are  numerous  possible  causes  of  transverse  angular  rates  of 
free-f light  rockets.  Some  are: 

1.  Exhaust  flow  interference,  including  "blow-by"  and  blast 
impingement  from  a  previously  launched  rocket  on  the  launcher. 

2.  Dynamic  imbalance  of  the  rocket. 

3.  Rocket  thrust  anomalies  due  to  mechanical  misalignment  of 
the  rocket  nozzle  and  unsteady  flow  within  the  rocket. 

4.  Elastic  deformations  of  the  rocket. 

5.  Elastic  deformations  of  the  launcher. 

6.  Guidance  mechanism  imperfections. 

Forces  caused  by  flow  interference  in  the  form  of  blow-by  may  be 
sufficient  to  produce  transverse  angular  rates  on  the  order  of  100  mrad/sec 
according  to  Booker.4  Dynamic  imbalance  can  cause  a  long  period  preces- 
sional  motion  of  the  unbalanced  rocket  which  appears  in  attitude  data  to  be 
a  secular  transverse  angular  rate  of  considerable  magnitude,  say  57  mrad/sec 
for  a  spin  rate  of  9  rps  and  a  1  mrad  angular  misalignment  of  the  principal 
axis  of  least  inertia  (PALI) .  Thrust  misalignment  of  a  mechanical  nature  may 
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produce  similar  rates.  As  mentioned,  elastic  vibration  of  a  spinning 
flexible  rocket  may  result  in  rates  on  the  order  of  100  mrad/sec,  while 
elastic  deformations  of  the  launcher  itself,  especially  the  launch  rail 
or  tube,  may  produce  similar  magnitude  rates. 

In  the  effort  reported  here,  elastic  deformations  of  the  rocket 
were  considered  of  principal  interest.  For  purposes  of  comparison  and  to 
match  optical  lever  data,  we  also  consider  thrust  misalignment.  We  note 
that  thrust  misalignment  may,  in  fact,  be  caused  by  a  static  elastic  de¬ 
formation  of  the  rocket.  Hence,  thrust  misalignment  and  elastic  deforma¬ 
tions  of  the  rocket  are  sometimes  related.  Also,  the  effects  of  dynamic 
imbalance,  which  could  be  caused  by  static  bending  of  the  rocket,  are 
much  like  those  of  thrust  misalignment.  If  the  thrust  and  spin  rate  of 
the  rocket  are  essentially  constant,  or  there  is  available  attitude  data 
over  only  a  short  time  period,  the  two  are  generally  indistinguishable. 

With  these  thoughts  as  background,  we  now  consider  transverse  vibrations 
as  "causes"  of  transverse  angular  rates. 

2.3  How  Transverse  Vibration  Causes  Transverse  Angular  Rates 

A  physical  explanation  for  the  way  in  which  transverse  vibration  of  a 
rocket  during  the  guidance  phase  results  in  a  transverse  angular  rate  of 
the  rocket  at  EOG  is  first  given  in  what  follows.  On  the  basis  of  this 
physical  explanation,  the  two-body  physical  and  mathematical  models  are 
used  to  obtain  a  simple  mathematical  expression  for  the  transverse  angular 
rate  which  results  from  transverse  vibration. 
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Physical  Explanation 

Assume  that  a  rocket  is  vibrating  transversely  prior  to  EOG.  The 
rocket  then  has  angular  momentum  due  to  motion  about  its  center  of  mass. 

This  angular  momentum  has  a  transverse  component  (see  Fig.  3)  that  is 
time  varying.  It  is  the  resultant  of  the  transverse  angular  momenta  of 
infinitely  many  segments  of  the  rocket  (considered  as  a  continuous  flexible 
body)  due  to  rotation  of  each  segment  about  its  center  of  mass  and  of  the 
angular  momenta  due  to  the  motion  of  the  segments'  centers  of  mass  relative 
to  the  center  of  mass  of  the  rocket. 


X 


hg  angular  momentum 
due  to  spin. 


b  angular  momentum 
1  about  C  due  to  transverse 
motion  ot  C . 


Fig.  3.  Angular  momentum  of  a  spinning  flexible  rocket. 

The  angular  momentum  of  each  segment  due  to  rotation  about  its  center 
of  mass  has  a  transverse  component  because  each  element  is,  in  general, 
rotating  transversely  to  the  line  segment  which  forms  a  portion  of  the 
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centroidal  axis  of  the  rocket.  It  also  normally  has  a  component  trans¬ 
verse  to  the  launch  axis  which  is  due  to  the  fact  that  angular  momentum 
resulting  from  spin  is  (because  the  rocket  is  bending  transversely)  not 
directed  along  that  axis. 

If  the  rocket  is  slender  and  is  not  bent  severely,  the  resultant  of 
the  angular  momenta  of  the  rocket  segments  about  their  respective  centers  of 
mass  is  much  less  than  that  due  to  the  motion  of  the  segments'  centers  of 
mass  relative  to  the  rocket's  center  of  mass.  The  latter  can  be  (as  will 
be  seen  later)  a  substantial  part  of  the  total  angular  momentum  of  the 
rocket. 

Until  EOG,  the  launcher  constrains  the  rocket  in  such  a  manner  that  the 
transverse  component  of  its  angular  momentum  due  to  transverse  (to  the  launch 
axis)  motion  changes  in  magnitude  and  direction.  Hence,  it  is  possible  that 
at  EOG  there  will  be  little  transverse  angular  momentum.  This  is,  however, 
not  very  probable.  Assuming,  then,  that  such  a  transverse  angular  momentum 
exists  at  EOG,  it  will  be  present  immediately  after  EOG  and  will  generally  be  at 
least  partially  transverse  to  the  instantaneous  principal  axis  of  least 
inertia  of  the  vibrating  rocket. 

Because  in  free-flight  the  rotational  and  vibrational  motions  of  the 
rocket  are  virtually  uncoupled  (i.e.,  it  behaves  somewhat  like  a  free-free 
beam) ,  if  the  principal  axis  of  least  inertia  of  the  rocket  is  rotating 
transversely  at  EOG,  it  will  continue  to  do  so  until  an  external  torque 
stops  it.  This  rotation  is  a  precession  and  is  periodic  if  no  external 
torques  are  present  and  there  is  no  energy  dissipation.  However,  the  period 
of  the  precession  is  very  long  for  slender  rockets,  since  it  is  approximately 
equal  to  2tt  1^/ (fi^I^),  where  I  is  the  transverse  moment  of  inertia  of  the 
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rocket,  and  is  its  spin  rate  and  1^  is  the  rocket's  axial  moment  of 
inertia.  Hence,  it  appears  to  be  constant  near  EOG. 

An  important  point  is  that  the  transverse  rate  due  to  transverse 
vibration  of  the  rocket  is  not  created  subsequent  to  EOG  as  a  rate  due  to 
blow-by  would  be.  It  is  the  result  of  energy  being  input  into  the  rocket 
for  a  period  of  time  before  EOG,  perhaps  from  the  time  of  ignition. 

Another  important  point  is  that  if  the  transverse  vibration  is  such  that  the 
PALI  of  the  vibrating  rocket  is  not  rotating  transversely  at  EOG,  then  there 
is  essentially  no  transverse  angular  rate  of  that  axis  due  to  the  vibration 
subsequent  to  EOG. 

Approximate  Transverse  Angular  Rate  Due  to  Vibration 

An  approximate  expression  for  the  transverse  angular  rate  at  EOG 
which  is  the  result  of  transverse  vibration,  is  derived  in  Appendix  C. 

There,  the  two-body  physical  model  and  the  constraints  used  in  the  two-body 
mathematical  model  are  used  to  determine  the  transverse  components  of  the 
angular  momentum  of  the  centroidal  principal  axes  of  the  rocket  at  EOG  in 

terms  of  model  physical  characteristics  and  the  pitch  and  yaw  rates  of  the 

•  • 

rocket's  nose  at  end  of  guidance,  0^  and  ¥  ,  respectively.  Note  that  in 


this  model,  and  are  due  to  bending  only. 

NO  NO  - 

Let  and  be  the  axial  and  transverse  moments  of  inertia, 
respectively,  of  Body  2;  a  =  m^m^/m, where  in.  is  the  mass  of  Body  j  and 
m  =  m^+m2  «  ri>  r2  and  ^2  tBe  axial  components  of  the  vectors  from 


P  to  C^,  Cx  to  R,  R  to  C2  and  to  Q,  respectively.  Then, 


“t  '  (6kl  ‘  W^NO  +  • 


(1) 


where  <5  =  tB2  +  0r2^1+r2^-Bl^IT~IA^  ’  kl  =  k2  “  d/(^2+1V  * 
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and  d  =  +  r^  +  +  r2‘  note  here  that  the  curves  on  page  37  of 

Ref.  2  are  not  all  correct,  since  incorrect  initial  conditions  were  used  to 
obtain  them.  The  correct  initial  conditions  are  given  in  Appendix  A. 

Equation  (1)  implies  that  is  directly  proportional  to  the  transverse 
angular  rate  of  the  rocket's  nose.  Regardless  of  the  transverse  angular 
rate  of  the  rocket's  nose,  u>t  will  be  zero  when  the  factor  6k^  -  k^/k^  is 
zero.  This  may  be  the  case  if  either  the  mass  and  inertia  characteristics 
of  the  two  bodies  are  altered  so  that  6  =  1/k^,  or  the  support  points  of 
the  rocket  are  changed.  The  latter  possibility  will  be  considered  shortly. 
First,  we  shall  see  what  size  transverse  angular  rate  may  be  expected. 

From  the  physical  characteristics  of  the  rocket  GEM  #7  (see  Table  1, 
infra)  the  following  values  needed  to  determine  the  constant  factor  in 
Eq.  (1)  can  be  obtained  (see  Table  2,  infra):  Bj  =  0.29  kg-m2 ,  B2=  11.86  kg-m2 , 
a  =  28.37  kg,  1^  =  0.406  kg-m2,  IT  =  99.7  kg-m2,  k^  =  1.22, and  k^  =  5.64. 

The  constant  factor  is  approximately  0.32.  Bending  oscillations  of  ampli¬ 
tudes  up  to  0.5  milliradians  have  been  observed  by  optical  lever  techniques. 
Assuming  such  an  amplitude  and  a  first  mode  bending  frequency  of  60  Hz,  we 
get  cot  =  60  mrad/sec,  a  substantial  transverse  angular  rate. 

If  we  replace  the  second  factor  in  Eq.  (1)  with  Au>v>  where  A  is  the 
amplitude  of  the  transverse  vibration  just  prior  to  EOG  and  where  U)v  is 
the  dominant  (probably  first-mode)  bending  frequency,  then  we  can  plot 
oj^/Auj^  as  a  function  of  ^ote  t^vat  ^  ma^  ne8ative?  indeed,  it 

is  in  the  above  example!)  The  curve  (see  Fig.  4)  obtained  is  virtually  a 
straight  line,  which  has  a  zero  crossing  at  ~  Thus, if 

£ 2  *  0.73  r2  had  been  used  in  the  above  example,  the  transverse  rate 
predicted  by  Eq.  (1)  would  have  been  approximately  zero. 


Fig.  4.  Transverse  angular  rate  as  a  function  of  forward 
support  location. 

Although  Eq.  (1)  does  not  include  the  effects  of  higher-order  bending 
modes,  it  still  yields  reasonably  good  estimates  of  the  transverse  angular 
rate  which  may  occur.  It  does  include  the  effects  of  support  placement  and 
the  results  presented  emphasize  the  importance  of  choosing  the  locations 
of  supports  judiciously. 

2.4  Results  from  Mathematical  Models 

To  illustrate  the  type  of  post-launch  attitude  motion  which  is  caused 
by  transverse  vibration  prior  to  EOG,  we  present  here  some  results  obtained 
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using  the  two-body  and  assumed  nodes  mathematical  models.  In  addition  to 
transverse  vibration,  the  effects  of  thrust  misalignment  on  a  flexible 
rocket  are  illustrated. 

Two-3ody  Model  Results 

Only  a  few  two-body  model  results  are  presented,  because  others  have 
been  given  in  Ref.  2. '  The  rocket  modeled  is  GEM  #7  (see  Table  1). 

Table  1.  Physical  Characteristics  of  GEM  ft 7 


Length  (m) 

3.353 

Radius  at  aft  end  (m) 

0.076 

Mass  (kg) 

114. 76 

IT  (kg-m2) 

99.7 

lx  C*g-®2) 

0.406 

Thrust  (Nt) 

46720 

Frequency  of  first  free- 
free  bending  mode  (Hz) 

=67 

The  data  used  in  the  two-body  model  is  given  in  Table  2.  It  should  be 
pointed  out  chat  chi3  data  is  such  that  the  center  of  mass  of  Che  undeforned 
rocket  model  is  located  at  the  hinge  point  R. 


Table  2.  Characteristics  of  Two-3ody  Model  of  GEM  -t 7. 


=L  (kg) 

50.466 

m7  (kg) 

64.286 

k  (Nt-m) 

1.49 1x10 5 

^  (m) 

0.9376 

-0.3354 

Ft  (NO 

46720 

r.  (m) 

X 

0.9376 

r2  (m) 

0.7393 

I,  (kg-m2) 

99.7 

Radius  (m) 

0.095 

Radius  (m) 

0.095 

IA  (kg-m2) 

0.406 

A1  (kg-m2) 

0.1631 

3,  (kg-m2) 

X 

0.1789 

(2^  (rad /sec) 

47 

A2  (kg-m2) 

10.996 

3?  (kg-m2) 

5.274 

Although  the  results  given  in  Ref.  2  are  not  exactly  correct  due  to  the  error 
in  initial  conditions  mentioned  previously,  the  error  had  very  little  affect 
on  the  results  given  there,  except  for  chose  on  page  37. 
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If  there  is  no  transverse  vibration  at  EOG,  but  there  is  thrust 
misalignment  defined  (see  Fig.  6)  by  o i  =  0.677  mrad  and  a.^  =  4.6  mrad, 
the  time  histories  shown  in  Fig.  7  are  obtained.  The  thrust  misalignment 
causes  only  a  small  amount  of  bending  so  that  the  curves  look  essentially 
like  those  obtained  using  a  rigid-body  model  for  the  rocket.  The  slope 

of  the  mean  value  of  0„  is  about  -80  mrad/sec. 

N 


Fig.  6.  Thrust  misalignment  angles. 

The  combined  effects  of  transverse  vibration  and  thrust  misalignment 
are  illustrated  by  the  curves  in  Fig.  8.  To  obtain  these  curves,  the  initial 
conditions  used  to  generate  the  results  shown  in  Fig.  5  and  the  thrust  mis¬ 
alignment  used  to  get  the  curves  in  Fig.  6  were  combined.  It  is  clear 
from  the  3^  curve  in  Fig.  7  that  the  transverse  rate  due  to  vibration  is 
more  significant  than  that  due  to  a  4  mrad  thrust  misalignment. 
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Assumed-Modes  Model  Results 

At  this  point  we  are  considering  the  effects  of  transverse  vibration 
regardless  of  the  cause.  Hence,  for  purposes  of  illustration,  a  bending 
oscillation  of  the  rocket  model  was  induced  while  it  was  "on  the  launcher" 
and  Qjj  and  4*  as  functions  of  time  were  generated  by  numerically  integrating 
the  assumed-mode  model  equations  of  motion.  The  data  used  is  given  in 
Table  3. 


Table  3.  Data  for  Assumed-Modes  Model  of  GEM  #7 


L  (m) 

3.353 

El  (Nt-m2) 

1.652  x  106 

a  (kg/m) 

(  26.93,  0  1  x  £  xM 
(  43.49,  x^j  <  x  <  L 

XM(m) 

1.875 

xQ<m)  '  2 

2.236 

I A  (kg-m2) 

0.406 

IT  (kg-m2) 

99.7 

Ft  (Nt) 

46720  (0.0  <  t  <  0.25) 

Ta  (Nt-m) 

284.2  (0.0  <  t  <  0.0647) 

End  of  guidance  occurred  (see  Fig.  9)  at  about  0.0647  seconds,  as  the 
nose  of  the  rocket  model  was  pitching  up.  The  transition  from  on- the- launcher 
motion  to  free-flight  is  characterized  by  (1)  a  change  in  the  slope  of  the 
mean  value  of  ©^  and  (2)  changes  in  the  amplitude  and  dominant  frequency 
of  the  bending  oscillations.  Such  changes  are  observable  in  optical  lever 
data.  The  natural  frequencies  of  the  first  bending  mode  of  the  model  for 
this  example  are  approximately  41  Hz  during  guidance  and  70  Hz  during  free 
flight. 
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TIME  (msec) 

Fig.  9.  Assumed-modes  results  -  effects  of  bending. 


The  effect  of  thrust  misalignment  is  qualitatively  the  same  when  the 

two-body  or  assumed-modes  model  is  used.  The  time  histories  shown  in  Fig. 

10  were  generated  by  using  a  and  a  values  stated  above  in  connection  with 

y  z 

the  two-body  model.  The  curves  shown  in  Fig.  10  differ  from  those  shown 
in  Fig.  7,  primarily  because  the  assumed-modes  model  of  the  rocket  had 
rotated  through  approximately  60  degrees  at  EOG.  Hence,  the  corresponding 
two-body  thrust  misalignment  angles  are  a  =  4.32  mrad  and 

ot^  -  -  1.72  mrad.  (The  roll  angle  in  the  two-body  mathematical  model  is 

_ A 


zero  at  EOG.)  When  the  proper  thrust  misalignment  angles  are  used  in  the 

two-body  model  and  the  pitch  and  yaw  angles  and  rates  of  the  rocket  model's 

nose  at  EOG  (as  taken  from  the  assumed-modes  model  results  with  vibration 

present)  are  used,  the  results  for  and  0„  as  functions  of  time  obtained 

N  N 

using  the  two  models  agree  very  well  (see  Fig.  11).  Since  the  computer 
codes  used  to  obtain  the  results  given  were  developed  at  different  points 
in  time  and  are  different  in  many  respects,  the  agreement  achieved  implies 
that  both  codes  are  correct. 


TIME  (msec) 

Fig.  10.  Assumed-modes  results  -  effects  of  thrust  misalignment. 
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The  rates  at  which  0„  and  f  are  changing  at  EOG  determine,  according 
N  N 

to  Eq.  (1),  the  transverse  angular  rate  of  the  rocket  after  EOG  which  is 
due  to  bending.  To  see  if  this  is  true  for  the  assumed-modes  model,  the 
guidance  length  was  changed  so  that  EOG  occurred  at  times  when  0^  was 
negative,  essentially  zero ,  and  positive.  The  results  for  0^(t)  are  shown 
in  Fig.  12.  They  are  in  agreement  with  the  two-body  model  results. 


TIME  (msec) 


Fig.  12.  Assumed-modes  results  -  effects  of  bending  with 
different  guidance  lengths. 

A  constant  factor  analogous  Co  that  in  Eq.  (1)  can  be  obtained  from 

assumed-modes  model  results  by  measuring  the  mean  slopes  of  the  0^(t)  an^ 

,(t)  curves  after  EOG  and  using  the  values  of  0„  and  V  at  EOG.  For  this 
N  N  N 

model  of  GEM  #7,  we  get  the  factor  0.33  which  is  in  close  agreement  with 
the  two-body  value  of  0.32. 
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As  a  final  illustrative  result,  we  offer  Fig.  13,  which  depicts  the 
pitch  and  yaw  angles  of  a  flexible  spinning  rocket  represented  using  the 
assumed-modes  model.  This  plot  of  0^  vs.  ¥  can  be  considered  as  "simu¬ 
lated"  optical  lever  data.  To  obtain  this  result,  we  used  the  data  in 
Table  3,  forced  0„  to  be  -1  mrad  and  'F  to  be  -120  mrad/sec  at  t=0,  by 
picking  the  proper  first  mode  generalized  coordinate  and  the  proper 
value  of  the  time  rate  of  change  p^  (see  Appendix  B  for  details) .  We  also 
set  a  *  -0.5  mrad  and  0.^  =  1  mrad.  Although  at  this  point  no  attempt  was 
made  to  match  any  actual  data,  the  "simulated"  data  is  in  some  ways  quali¬ 
tatively  similar  to  that  actually  obtained  during  the  flight  of  GEM  #8 
(see  Fig.  14).  However,  the  complexity  of  the  actual  optical  lever  data 
shows  bending  modes  other  than  the  first  were  excited. 


Fig.  13.  Assumed-modes  results  -  "simulated"  optical  lever  data. 


SECTION  3.  ANALYSIS  OF  OPTICAL  LEVER  DATA 


3.1  Introduction 

The  optical  lever  is  the  most  accurate,  commonly-used  means  of 
obtaining  attitude  time  histories  of  free-flight  rockets.  Optical  lever 
data  is  usually  acquired  during  the  guidance  phase  and  a  short  portion  of 
the  free-flight  phase.  Through  appropriate  processing,  the  raw  data  is 
converted  into  pitch  and  yaw  angles  of  the  rocket's  nose  as  functions  of 
time.  If  the  rocket  is  not  vibrating  transversely  during  the  data 
acquisition  time,  the  pitch  and  yaw  angles  of  the  rocket's  nose  are,  of 
course,  those  of  the  rocket.  However,  if  it  is  vibrating  transversely, 
this  is  clearly  not  so.  In  fact,  if  the  rocket  is  vibrating  transversely, 
there  arises  the  question  of  what  coordinate  frame  should  be  used  in 
specifying  the  "attitude  of  the  rocket." 

The  coordinate  frames  most  commonly  used  in  specifying  the  attitude 
of  an  unconstrained  flexible  body  are:  (1)  a  coordinate  frame  rigidly 
fixed  in  an  essentially  rigid  portion  of  the  body,  (2)  principal  axes  of 
the  deformed  body, and  (3)  so-called  "mean  axes"  which  are  such  that  the 
angular  momentum  with  respect  to  the  mean  axes  due  to  vibration  is  zero 
at  all  times.  As  mentioned  previously,  principal  axes  of  the  deformed 
rocket  are  used  in  the  assumed-modes  mathematical  model  during  the  free- 
flight  phase. 

If  the  transverse  vibrations  are  assumed  to  have  no  effect  on  the 
rocket's  gross  attitude  motion  subsequent  to  EOG,  the  data  can  be  analyzed 
by  assuming  that  the  variation  in  the  attitude  of  the  rocket  is  due  to: 
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(1)  an  initial  transverse  angular  velocity  due  to  unspecified  factors, 

(2)  the  action  of  a  body-fixed  torque  due  to  thrust  misalignment, 

(3)  dynamic  imbalance  of  the  rocket,  and  (4)  space-fixed  torques  arising 
from,  for  example,  blow-by  and  gravity  (if  a  tip-off  launcher  is  used). 

If  the  dynamic  imbalance  of  a  rocket  which  is  spinning  at  a  constant  rate 
is  also  constant,  it  produces  variations  in  rocket  attitude  which  are  of 
the  same  form  as  those  produced  by  an  "equivalent"  constant  body-fixed 
torque.  Thus,  if  a  constant  spin  rate  model  is  used  in  the  analysis 
process,  constant  dynamic  imbalance  and  constant  thrust  misalignment 
effects  are  not  distinguishable.  Furthermore,  since  the  torque  due  to 
blow-by  acts  approximately  impulsively,  it  is  difficult  to  distinguish  its 
effects  from  those  of  the  otherwise  present  initial  angular  velocity.  For 
these  reasons,  in  the  analysis  presented  here,  only  the  effects  of  initial 
transverse  angular  velocity  due  to  unspecified  factors  and  thrust  mis¬ 
alignment  are  considered  when  the  effects  of  transverse  vibrations  are 
ignored.  When  the  effects  of  transverse  vibration  are  not  ignored,  the 
initial  (at  EOG)  angular  velocity  due  to  it  and  possibly  an  additional 
initial  angular  velocity  due  to  other  factors  (launcher  vibration,  for 
example)  as  well  as  thrust  misalignment  is  considered.  The  effects  of 
blow-by  are  not  considered  because  pressure  data4  taken  during  the  guidance 
phase  of  the  rocket  of  interest  indicated  blow-by  was  not  significant. 

3.2  Transverse  Vibration  Effects  Ignored 

If  the  rocket  is  modeled  as  a  rigid  body,  which  is  spinning  at  a 
constant  rate  and  is  acted  upon  by  a  constant  body-fixed  torque  due  to 
a  constant  magnitude  thrust,  F^,,  that  is  not  aligned  with  the  longitudinal 
axis  (x-axis)  of  the  rocket,  Euler's  equations  and  Poisson's  kinematic 
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equations  can  be  solved  to  obtain  closed-form  solutions  for  the  yaw, 

pitch  and  roll  angles  of  its  centroidal  principal  reference  frame. 

Furthermore,  if  the  rocket  is  slender,  so  that  I  «  I_,  and  the  pitch 

A  I 

and  yaw  angles  are  small,  the  solutions  for  0^  and  V  can  be  well  approx¬ 
imated  for  a  short  period  of  time  (say  one-half  a  second)  by  the  expres¬ 
sions  given  on  Fig.  15.  In  Fig.  15,  a  and  a  are  "pitch"  and  "yaw"  angles 

y  z 

of  the  thrust  vector  with  respect  to  the  rocket's  x-axis.  Also,  and 
are  the  initial  transverse  angular  velocity  components.  The  curves  in 
Fig.  16  are  typical  ones.  Certain  geometric  characteristics  of  the  curves 


are  defined  in  Fig.  15.  These  are  expressed  analytically  in  terms  of 

(2  ,  fl  (2  I  ,  I  ,  a  ,  a  ,  F„  and  &  ,  the  distance  from  the  intersection 
x  yO  zO  T  A  y  z  T  C 

of  F^  with  the  x-axis  and  the  rocket's  center  of  mass.  By  labeling  plots 

of  0„  and  ’V  versus  time  obtained  from  optical  lever  data  as  shown  in 
N  N 

Fig.  15,  if  Q.  ,  F  ,  i  ,  I  and  I  are  known,  ct  ,  a  ,  0  „  and  Q  _  can  be 
xTcT  A  yzyO  zO 

obtained  graphically. 

To  give  an  example  of  this  "graphical  analysis,"  we  use  optical 


lever  based  data  for  GEM  #8.  The  physical  characteristics  of  GEM  #8  are 


given  in  Table  4.  From  the  appropriately  labeled  curves,  we  get: 

8  -  -18  mrad/sec,  A  -  0.945  mrad,  8  -  20  mrad/sec,  A  -  2.41  mrad, 

y  y  z  z 

and  A  =  2.6  mrad. 

For  GEM  #8,  IT  -  99.5  kg-m2 ,  I  -  0.421  kg-m2 ,  l  -  1.874  m, 

x  A  C 

4* 

J2  =*  30  rad/sec  and  FT  =  40940  Nt.  Hence,  K  -  0.8642.  It  follows  from 

the  definitions  of  8  ,  A  ,  8  and  A  that  fi  _  =  -90  mrad/sec, 

y  y  z  z  yO 

(2  .  -  48.5  mrad/sec,  a  =  1.104  mrad  and  a  -  2.79  mrad. 

zO  y  z 


This  is  an  average  value  of  fljj,  since  GEM  #8  was  spun  via  an  eroding  spin 
turbine.  It  is  actually  lower  than  a  simple  average,  because  the  nutational 
frequency  does  not  instantaneously  become  (I  -  / I  . 


E 


<£=£»  *Aco*(£y-y)  -4 


TIME 


©=/Syt-A*ln(fl,t-y)-4y 

y=tan'(av/aI) 


^8*  Pitch  and  yaw  angles  when  thrust  misalignment  is  present 


Table  4.  Physical  Characteristics  of  GEM  #8 


Length  (m) 

3.353 

Radius  at  aft  end 

(m) 

0.076 

Mass  (kg) 

114.5 

IT  (kg-m2) 

99.5 

IA  (kg-m2) 

0.421 

Thrust  (Nt) 

40940 

Frequency  of  first  free- 
free  bending  mode 

-  67 

Spin  rate  at  EOG 

(rad/sec) 

-  30 

Distance  from  aft 
initial  center  of 

end  to 
mass  (m) 

1.874 

0 


50 


200 


250 


—i — I - r— 

100  150 

TIME  (msec) 

Fig.  16.  Comparison  of  theoretical  results  and  data  -  bending 
effects  ignored. 
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If  the  above  values  of  ft  ft  n,  a  and  a  are  substituted  into  the 

yO  zO  y  z 

expressions  for  @N  and  4^  shown  in  Fig.  15,  the  solid  curves  in  Fig.  16 
are  obtained.  The  dashed  curves  in  Fig.  16  are  those  derived  from  optical 
lever  data.  Very  close  agreement  between  the  theoretical  and  experi¬ 
mental  curves  is  apparent  as  far  as  the  "mean"  values  of  0^  and  4*  are 
concerned,  even  though  the  model  is  one  with  constant  spin  rate. 

The  magnitudes  of  the  values  of  ft^  and  ft^  required  to  obtain  the 
"match"  illustrated  in  Fig.  16  are  very  large.  Since  blow-by  did  not 
appear  to  be  a  significant  causative  factor  in  the  case  of  GEM  #8,  the 
cause  of  these  large  magnitude  initial  angular  velocity  components  must 
be  found  elsewhere.  In  the  next  subsection,  we  show  that  it  is  very 
probable  that  one  cause  was  transverse  vibration  during  guidance. 

3.3  Transverse  Vibration  Effects  Included 


To  include  the  effects  of  transverse  vibration  in  the  analysis  of 


optical  lever  data,  we  use  the  two-body  mathematical  model.  Two  approaches 

are  taken.  In  the  first  we  :  (1)  graphically  obtain  0  ^,  4>  ,  0  and  4'.,. 

nu  NO  NO  NO 

estimates  from  the  curves  derived  from  optical  lever  data,  (2)  use  0^ 

and  4'Nq  in  Eqs.  (C-7)  (see  Appendix  C)  to  get  and  due  to  the 

vibration  and  (3)  add  to  and  oj^q  initial  angular  velocity  components 

due  to  "other"  factors  to  obtain  the  ftyQ  and  found  in  subsection  3.2. 

The  second  approach  is  to  fit  the  optical  lever  derived  curves  with  the 

two-body  model  solution  in  a  least-squares  sense. 

Taking  the  first  approach,  we  first  find  0^Q  =  0.25  mrad,  4*  =  0.0 

mrad,  0.Trt  -  -220  mrad/sec  and  4'  =  140  mrad/sec.  For  the  GEM  #8  two- 
NO  NO 

body  model  B^~  0.194  kg-m2,  B^  -  5.088  kg-mz,  m^  =  50.51  kg  and 

~  64.35  kg.  Hence,  a  =  28.22  kg  and  the  constant  factor  (6k^-k^/k^)  ■ 
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0.2733.  From  Eqs.  (C-7),  c^q  -  "60  mrad/sec  and  0)^  =  38  mrad/sec  due  to 
transverse  vibration.  The  total  transverse  rate  due  to  bending  is, 
therefore,  approximately  70  mrad/sec. 

Care  must  be  exercised  to  maintain  consistency  when  using  the  graph¬ 
ical  method.  Note  that  A  =  /A*  +  A2  ;  so, the  values  of  A  and  A  and  A 

y  z  y  z 

must  be  consistent.  Also,  the  equation.Y  =  tan-1 (A  /A  ),  must  be  approx- 

y  z 

imately  satisfied  by  the  measured  value  of  y.  Attaining  such  consistency 
generally  requires  several  attempts  at  drawing  the  enclosing  lines  and 
centerlines  of  the  oscillations. 

It  is  obviously  no  small  problem  to  determine  exactly  what  f  , 


0^  and  'Fjjq  are  from  the  optical  lever  data.  Again,  the  theoretical 
results  must  be  consistent  in  that  the  amplitudes  of  the  bending  oscillations 
in  the  theoretical  results  must  match  well  with  those  in  the  optical  lever 
data.  The  results  shown  in  Fig.  17  (see  page  33)  do,  in  fact,  agree  well. 

The  second  approach  involves  considerably  more  labor  than  the  graph¬ 
ical  approach, unless  the  optical  lever  data  is  available  in  a  form  suitable 
for  input  into  a  digital  computer.  It  is  a  standard  least-squares  procedure 
in  which  influence  functions  obtained  from  the  two-body  model  computer  code 
are  used.  The  influence  functions  are  found  by  setting  all  but  one  of 
the  initial  conditions  and  thrust  misalignment  angles  successively  to  zero. 
The  nonzero  quantity  is  set  equal  to  the  same  constant  in  each  case.  We 
used  the  values  of  1  mrad  for  angles  and  1  mrad/sec  for  rates  to  obtain  the 
influence  functions  used  to  get  the  results  which  follow  in  this  section. 

To  put  the  method  in  mathematical  terms,  we  let  ^  be  a  nxl  matrix 

of  angles  extracted  from  optical  lever  data.  We  also  let  x  *  (0,T/. 

-  NO  NO  NO 

•  T 

a  a  )  and  A  denote  a  nx6  matrix  of  influence  coefficients,  i.e., 

NO  y  z  ■  - 


jit 
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values  of  the  influence  functions  evaluated  at  times  after  EOG  corre¬ 
sponding  to  the  data  points.  Then,  if  the  two-body  model  were  assumed 
to  give  exact  results,  we  would  have 

l  =  A  x  .  (2) 

Since  the  two-body  model  cannot  be  reasonably  expected  to  give  exact 
results  and  since  we  would  like  to  use  more  than  six  pieces  of  data, 
we  let  the  error  involved  be 

<  =  y  -  A  x  .  (3) 

We  then  choose,  as  our  "answer"  for  x,  the  vector  (6x1  matrix)  x, 
such  that  £  €  is  minimized  when  x  is  substituted  for  x  in  Eq.  (3). 

The  vector  x  is  given  by  the  equation, 

x  =  (AT  4)-1  AT  1  .  (4) 

By  using  25  values  of  4*  ,  22  values  of  0  and  the  appropriate 

N  N 

influence  coefficients  from  the  two-body  model  computer  code  (in  which 

is  set  equal  to  30  rad/sec)  in  Eq.  (4),  the  following  values  may  be 

obtained:  0„-  =  0.8722  mrad,  0„-  *  -387  mrad/sec,  ?  =  0.1947  mrad, 

NU  NO  NO 

A 

=  128.5  mrad/sec,  a  *  0.685  mrad,  a  =  3.214  mrad.  These  estimates 
NO  y  z 

are  in  reasonably  close  agreement  with  those  obtained  using  the  graphical 
procedure.  One  difference  is  that  in  obtaining  the  former  we  ha>/e  not 
considered  causes  of  transverse  angular  rates  at  EOG  other  than  trans¬ 
verse  bending.  It  turns  out  that  if  one  attempts  to  solve  for  transverse 

T 

angular  velocity  components  due  to  "other"  factors,  the  matrix  A  A  cannot 
be  inverted. 

When  the  estimates  given  immediately  above  are  substituted  into  the 
two-body  model  computer  code,  the  results  shown  as  solid  lines  in  Fig.  18 
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are  obtained.  Good  agreement  between  the  theoretical  results  and  experi¬ 
mental  data  is  apparent.  However,  the  bending  oscillations  in  4*  are 
too  small  and  those  in  0N  are  too  large.  This  is  no  doubt  due  to  the 
presence  of  transverse  angular  velocity  components  at  EOG  due  to  factors 
other  than  transverse  bending. 

If  we  "adjust"  the  estimates  to  account  for  the  transverse  angular 
velocity  components  due  to  "other"  factors,  we  get  the  solid  curves  in 

A 

Fig.  19.  The  "adjusted"  values  are  =  0.8722  mrad,  0^Q  =  -208  mrad/ 
sec,  Qjq  =  0.1947  mrad  and  =  166  mrad/ sec.  In  addition,  transverse 

angular  rates  of  -49  mrad/sec  and  10  mrad/sec  in  pitch  and  yaw,  respec¬ 
tively,  are  required.  The  "fit"  is  very  good.  It  follows  that  the  least- 
squares  method  combined  with  a  little  "engineering  judgment,"  perhaps 
influenced  by  graphical  analysis  results  will  provide  useful  information 
concerning  the  effects  of  transverse  vibrations  and  "other  factors"  which 
produce  transverse  angular  rates  at  EOG. 
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SECTION  4.  PROBABLE  CAUSES  OF  TRANSVERSE  VIBRATION 


4.1  Introduction 

It  is  clear  from  the  results  given  in  the  previous  sections  that  if 
transverse  vibration  of  a  free-flight  rocket  is  occurring  at  EOG  there 
will  generally  be  an  associated  transverse  angular  rate.  The  exact  cause, 
or  causes,  of  the  bending  evidenced  in  GEM  series  flight  test  data  has 
not  been  determined  during  the  current  research  effort.  However,  several 
"probable"  causes  have  been  identified  and  studied  by  the  authors  and/or 
by  U.S.  Army  Missile  Laboratory  personnel  or  U.S.  Army  contractors. 

Soon  after  the  search  for  causes  of  transverse  bending  began,  it  was 
discovered  that  some  of  the  GEM  series  rockets  were  supported  during 
portions  of  their  guidance  phases  at  three  longitudinal  points  rather  than 
two  points.  Hence,  it  is  probable  that  this  "intermittent"  third  support 
caused  bending  of  the  rocket . 

A  second  cause  of  bending  was  suggested  by  Golden.5  As  the  rocket 
begins  to  spin  within  the  tube,  transverse  frictional  forces  act  at  support 
points.  These  forces  may  alternate  in  direction  and  hence  cause  bending. 

A  third  cause  of  bending  has  been  studied  briefly  by  Smith.6  If  the 
thickness  of  the  rocket  motor  case  is  not  uniform  circumferentially,  the 
pressurization  of  the  motor  case  during  ignition  of  the  rocket  motor,  and 
the  axial  forces  on  the  rocket  as  it  travels  on  the  launcher  will  generally 
cause  the  rocket  motor  case  to  bend.  Since  the  axial  forces  which  act 
during  detent  and  guidance  are  not  constant,  it  is  reasonable  to  believe 
that  transverse  vibration  of  the  rocket  with  an  eccentric  motor  case  will 
occur. 
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In  the  following  subsections,  we  treat  the  first  two  of  the  causes 
discussed  above.  We  first  consider  the  effects  of  "multiple"  (more  than 
two)  supports. 

4.2  Multiple  Supports 

The  computer  code  which  contains  the  assumed-modes  mathematical  model 
was  modified  so  that  three  points  of  support  could  be  modeled.  This  was 
done  by  forcing,  mathematically,  a  point  x^  on  the  centerline  of  the  rocket 
to  follow  a  prescribed  path  for  a  prescribed  amount  of  the  total  guidance 
length.  The  path  chosen  was  a  ramp.  That  is,  the  rocket  travels  within 
the  launch  tube  for  a  certain  distance.  Then,  the  point  x^  on  the  center- 
line  of  the  rocket  is  forced  to  follow  a  straight-line  path  which  is  not 
collinear  with  the  launch  axis,  i.e.,  the  ramp.  After  the  rocket  travels  a 
certain  additional  distance  within  the  launch  tube,  the  "third  support" 
ceases  to  support  the  rocket;  however,  the  other  two  supports  are  still 
present  until  EOG.  Hence,  the  rocket  is  forced  to  bend  and  is  then  released 
(partially) .  The  resulting  transverse  vibration  of  the  rocket  model  exists 
until- and  subsequent  to  EOG. 

Only  typical  results  are  presented  here  to  illustrate  the  possible 
effects  of  more  than  two  longitudinal  support  locations.  The  curves  in 
Fig.  20  were  obtained  using  the  data  given  in  Table  3.  The  sequence  of 
events  is  as  follows.  The  point  x^  corresponding  to  the  head  end  of  the 
motor  case  of  GEM  #7  is  forced  to  move  vertically  down  0.0254  cm  during 
0.1524  m  of  travel  in  the  x^-direction.  The  "third  support"  which  caused 
the  displacement  of  x^  is  then  removed.  The  induced  deflection  of  the  rocket 
model's  nose  is  about  0.8  mrad  when  the  third  support  is  removed.  The 
rocket  model  then  begins  to  vibrate  and  a  transverse  rate  of  about  61 
mrad/sec  is  present  at  EOG.  Hence,  a  small  displacement  of  Xj^  during 
guidance  is  sufficient  to  produce  a  large  transverse  rate  at  EOG. 


PITCH  ANGLE  (mrad)  »  YAW  ANCLE  (inrad) 
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4.3  Flexible  Supports  and  Friction 

The  two-body  model  of  a  spinning  flexible  rocket  was  extended  by 
constraining  the  two  bodies  elastically  during  the  guidance  phase  and 
adding  a  model  for  the  friction  between  the  rocket  and  the  flexible 
supporting  tube.  The  extended  model  is  described  in  Appendix  D.  A 
computer  code  was  written  to  solve  the  equations  of  motion  numerically. 
Numerical  solution  is  necessary  because  of  the  time-varying,  nonlinear 
character  of  the  equations  of  motion.  By  using  the  computer  code,  one  may 
study  the  effects  of :  (1)  flexibility  of  the  supports  (tube  and  shoe  flex¬ 

ibility),  (2)  friction  between  the  shoes  and  the  tube  wall,  (3)  clearance, 
or  "dead  space,"  between  the  shoes  and  wall,  (4)  forward  support  location, 

(5)  spin  torque  magnitude  and  (6)  tipoff.  Only  the  first  four  of  the  above 
will  be  discussed  in  what  follows. 

Flexibility  of  the  Supports 

If  frictionless  motion  is  assumed,  the  time  histories  0,T  and  V  shown  in 

N  N 

Fig.  21  are  obtained  using  the  physical  parameters  for  GEM  #8  and  the 
support  stiffness  Kp  =  =  8.76  x  107  Nt/m.  There  is  an  initial  deflection 

of  the  rocket's  nose  due  to  gravity,  but  until  EOG  little  change  in  this 
deflection  occurs.  At  EOG,  a  mean  transverse  rate'  of  only  0.67  mrad/sec 
is  present.  Variations  in  the  support  stiffness  produce  essentially  no 
difference  in  the  mean  transverse  rate  at  EOG. 

Friction 

When  non-zero  coefficients  of  friction,  Up  and  Uq,  are  introduced  into 
the  equations  of  motion,  the  rocket  generally  begins  to  "whirl"  as  the  spin 
rate  increases  and  the  frictional  forces  cause  transverse  vibration.  This 


Here,  "mean  transverse  rate"  is  the  square  root  of  the  sum.of  the. squares 

of  the  mean  values  (bending  oscillations  not  included)  of  0  and  tL. 

N  N 


motion  is  unstable;  i.e.,  the  bending  oscillations  increase  in  amplitude 
as  time  progresses.  Fig.  22  shows  the  time  histories  of  0^  and  ¥ 
obtained  using  the  above  stated  values  of  parameters  and  friction  coef- 
icients  dp  =  y^  =  0.05.  The  mean  transverse  rate  at  EOG  is  greater  than 
with  no  friction, but  is  still  less  than  0.8  mrad/sec. 

The  mean  transverse  rate  is  much  larger  for  the  case  of  yp  =  y^  =  0.1 
(see  Fig.  23).  Its  value  is  about  240  mrad/sec.  Larger  values  of  dp  and 
y^  thus  produce  more  bending  and  larger  mean  transverse  angular  rates.  In 
fact,  hp  and  y^  =  0.2  produce  divergent  motion  to  the  extent  that  the  small 
angle  assumption  is  violated  for  ©2  and  0^  after  only  0.05  seconds. 

When  friction  is  present,  the  stiffnesses  of  the  supports  become 
much  more  important.  For  Kp  =  =  4.38  x  107  Nt/m  and  yp  =  y^  =  0.05,  the 

mean  transverse  rate,  as  stated  above,  is  only  about  0.8  mrad/sec;  whereas, 
for  Kp  =  Kq  =  1.75  x  108  Nt/m  and  the  same  yp  and  y^  values,  it  is  approxi¬ 
mately  2  mrad/sec.  This  larger  rate  is  probably  generated  because  the 
"violence"  with  which  the  rocket  is  forced  to  and  fro  in  the  tube  when  the 
tube  is  very  stiff  causes  more  bending. 

Clearance 

If  the  tube  is  large  enough  that  all  the  shoes  (four  assumed)  at  each 
support  location  are  not  all  in  contact  with  it  initially,  there  will  be 
periods  of  time  during  guidance  in  which  the  rocket  is  not  supported  at 
one,  or  both,  of  the  support  locations.  The  distance  through  which  a  point 
on  a  centerline  of  the  rocket, at  a  support  point  and  originally  on  the 


centerline  of  the  tube, must  move  for  a  shoe  to  contact  the  tube  wall  is 


denoted  by  e.  In  the  work  we  have  done,  e  has  been  assumed  to  be  the  same 
at  the  fore  and  aft  support  points. 
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We  have  found  that  a  nonzero  e,  say  0.31  mm,  with  no  friction,  will 
result  in  a  slightly  greater  mean  transverse  rate  at  EOG  than  would  other¬ 
wise  be  present. 

An  example  of  the  attitude  motion  of  the  rocket's  nose  that  occurs 
when  Kp  =  =  8.76  x  107  Nt/m,  Up  and  =  0.05  and  £  =  0.31  mm  is  shown 

in  Fig.  24.  The  irregularity  of  the  time  histories  of  and  0^  during  the 
guidance  phase  is  due  to  the  clearance.  During  certain  intervals  of  time, 
the  rocket  is  actually  in  free-flight  within  the  tube. 

Forward  Support  Location 

If  the  forward  support  is  moved  from  its  nominal  =  ~0.335  m) 
location  on  the  rocket,  the  equilibrium  pitch  angle  of  the  nose  changes. 

For  the  two-body  model,  if  the  support  is  moved  so  that  =  0.54  m  (the 
value  found  in  Section  2  which  makes  the  constant  factor  of  Eq.  (1)  zero,  then 
the  equilibrium  pitch  angle  of  the  rocket's  nose  is  about  0.24  mrad.  For 
£2  =  0.54  m,  Kp  =  Kq  =  8.76  x  107  Nt/m,  =  0.05,  £  =  0  mm,  and 

all  other  values  of  parameters  the  nominal  ones,  we  get  the  result  shown 
in  Fig.  25.  The  mean  transverse  rate  is  not  zero,  but  is  approximately 
12  mrad/sec.  It  is  not  zero  because  the  support  points  P  and Q are  moving 
at  EOG  in  such  a  way  that  there  is  a  transverse  angular  rate  of  the  rocket's 
principal  axes  in  addition  to  the  rate  produced  by  transverse  vibrations. 
Hence,  the  suggested  method  for  reducing  the  adverse  effects  of  transverse 
vibration,  i.e.,  placing  the  supports  properly,  does  not  necessarily  pro¬ 
duce  zero  transverse  angular  rate  when  other  factors,  such  as  tube  flex¬ 
ibility,  are  present. 

Combined  Effects 

As  a  final  example  of  the  effects  of  friction  clearance  and  support 
location,  we  offer  Fig.  26.  The  results  shown  in  Fig.  26  were  obtained 
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using  ^2  =  0.54  m,  =  8.76  x  107  Nt/m,  =  0.05,  £  =  0.31  mm 

and  the  thrust  misalignment  angles  =  az  =  1  mrad.  In  this  case,  the 
transverse  rate  is  less  than  that  which  would  be  present  if  only  thrust 
misalignment  existed. 

The  qualitative  similarity  between  the  pitch  and  yaw  time  histories 
shown  in  Fig.  26  during  the  guidance  phase  and  the  optical  lever  results 
shown  in  Figs.  16,  17,  18  or  19  should  be  noted.  This  is  encouraging,  as 
it  indicates  that  a  mathematical  model  which  includes  the  flexibilities 
of  the  rocket  and  its  supporting  structure,  clearance,  friction,  thrust 
misalignment,  dynamic  imbalance  and  possibly  only  a  few  more  characteris¬ 
tics  could  be  used  successfully  in  analyses  of  free-flight  rocket  motion 
during  both  the  guidance  and  free-flight  phases. 

Comments  on  Angular  Momentum 

That  the  transverse  component  of  angular  momentum  may  represent  a  sub¬ 
stantial  portion  of  H,  the  total  angular  momentum  of  a  flexible  rocket  about 
its  center  of  mass  at  EOG  was  verified  by  using  the  extended  two-body  model. 
The  Body  1-fixed  (Cxyz  frame)  and  space-fixed  components  of  H  were  computed 
along  with  the  rocket's  state.  The  transverse  component  of  H  at  EOG  was 
found,  for  example,  to  be  0.633  | H j  for  the  combined  effects  case. 

The  space-fixed  angular  momentum  components  also  served  as  checks  on 
the  numerical  integration  procedure.  If  there  are  no  external  torques 
subsequent  to  EOG,  they  should  remain  constant.  Constancy  to  four  signif¬ 
icant  figures  was,  in  fact,  achieved. 


SECTION  5.  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Conclusions 


On  the  basis  of  results  of  this  investigation  we  have  reached  the 
following  conclusions: 

1.  The  two-body  model  and  assumed-modes  model  computer  codes 
are  correct  as  far  as  programming  is  concerned. 

2.  The  assumed-modes  computer  code  allows  one  to  model  a 
flexible  free-flight  rocket  in  more  detail  than  is  possible 
with  the  two-body  computer  code.  The  changes  in  amplitudes 
and  frequencies  of  transverse  vibrations  obtained  using  the 
assumed-modes  code  are  qualitatively  the  same  as  those  which 
are  often  observed  in  optical  lever  data. 

3.  The  graphical  method  presented  is  useful  in  getting  order 
of  magnitude  information. 

4.  The  two-body  model  is  a  valid  representation  of  a  spinning, 
flexible  free-flight  rocket  to  the  extent  that  it  can  be 
used  to  analyze  optical  lever  data  which  contain  evidence 
of  transverse  vibration. 

5.  The  methods  for  analyzing  optical  lever  data  which  are 
given  in  the  body  of  this  report  provide  useful  information 
concerning  the  effects  of  transverse  vibrations.  They  also 
can  be  used  to  detect  the  presence  of  "other"  factors  which 
cause  transverse  angular  rates  of  free-flight  rockets  at 
EOG. 

6.  Transverse  vibration  of  a  spinning  flexible  free-flight 
rocket  during  the  guidance  phase  can  result  in  transverse 
angular  rates  of  the  rocket  as  a  whole  at  EOG.  These  rates 
are  of  the  same  magnitude  as  those  caused  by  several  milli- 
radians  angular  thrust  misalignment,  generally  being  on  the 
order  of  10  to  100  mrad/sec. 

7.  The  transverse  angular  rate  due  to  bending  vibrations  can 
be  reduced  by  placement  of  the  rocket  supports  in  such  a 
way  that  rotational  motion  of  the  time-varying  centroidal 
principal  axes  of  the  rocket  is  reduced. 

8.  A  probable  cause  of  transverse  bending  of  some  of  the  GEM 
series  rockets  was  the  use  of  an  "intermittent"  third  support. 
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9.  Coulomb  friction  acting  axially  and  circumferentially  on 
a  spinning,  flexible  rocket  at  its  points,  or  areas,  of 
support  may  cause  considerable  transverse  bending,  the 
magnitude  of  which  depends  on  the  coefficient  of  friction, 
clearance  and  other  factors.  Friction  is  another  probable 
cause  of  bending  vibrations  of  the  GEM  series  rockets. 

10.  If  the  points  of  support  of  a  spinning,  flexible  rocket  are 
not  constrained  to  translate  in  a  fixed  direction  during 
guidance,  then,  although  the  transverse  rate  due  to  bending 
vibrations  may  be  reduced  by  choosing  the  support  points 
properly,  the  transverse  angular  rate  of  the  rocket's 
principal  axes  may  still  be  large  due  to  the  transverse 
motions  of  the  support  points . 

11.  Additional  work  is  needed  to  provide  a  complete  understand¬ 
ing  of  the  causes  and  effects  of  transverse  vibration  of 
spinning,  flexible  rockets  supported  and  guided  by  flexible 
launchers.  The  results  presented  herein  should  provide  a 
basis  for  such  work. 

Conclusions  1  and  2  are  based  on  the  results  of  comparisons  of  the 
two-body  and  assumed  modes  computer  codes.  Conclusion  3  is  an  obvious 
result  of  the  successful  use  of  the  graphical  method  in  analyzing  GEM  #8 
data  in  Section  3.  Conclusion  4  is  reached  in  light  of  the  success 
achieved  in  Refs.  2  and  3  and  herein  in  "fitting"  optical  lever  data  with 
the  two-body  results.  There  is  still  some  question  about  uniqueness  of 
the  parameter  set  needed  to  "fit"  the  data. 

The  analysis  presented  in  Section  3  justifies  Conclusion  5.  From 
results  given  in  Sections  2  and  3,  it  is  apparent  that  transverse  vibra¬ 
tion  during  guidance  is  a  significant  cause  of  transverse  angular  rates 
subsequent  to  EOG,  as  concluded  in  6,  above.  Results  given  herein  also 
indicate  that  if  the  rotational  motion  of  the  time-varying  centroidal 
principal  axes  of  a  flexible  rocket  is  due  only  to  bending  of  the  rocket, 
such  rotation  can  be  reduced  (Conclusion  7).  Conclusions  8,  9  and  10 
follow  from  the  results  of  Section  4. 


Although  it  is  felt  that  considerable  progress  has  been  made  since 
the  work  reported  in  Ref.  2  was  begun,  the  problem  of  determining  the 
motion  of  a  general,  flexible  free-flight  rocket  while  it  is  being  sup¬ 
ported  and  guided  by  a  general  flexible  launcher  has  obviously  not  been 
solved.  A  complete  solution  is  apparently  not  attainable.  However,  as 
Conclusion  11  indicates,  much  more  can  and  needs  to  be  done.  We  consider 
this  next. 

5 . 2  Recommendations 

Several  problem  areas  should  be  more  carefully  explored  in  the  near 


future.  These  are  considered  in  the  following  recommendations. 


1.  Additional  work  is  recommended  with  the  goal  of  developing  a 
more  accurate  and  efficient  method  for  correlating  theoret¬ 
ical  and  experimental  results  for  spinning , flexible  free- 
flight  rocket  motion.  More  general  mathematical  models  for 
both  the  launcher  and  rocket  should  be  developed  so  that  the 
simultaneous  effects  of  launcher  motion,  launcher  flexibility 
and  aerodynamic  effects  (including  blow-by)  can  be  studied. 
Such  models  coupled  with  a  statistical  parameter  estimation 
algorithm  would  be  very  useful  in  analyzing  test  results. 


2.  Although  it  is  now  agreed  by  many  that  transverse  vibration 
of  a  free-flight  rocket  during  guidance  does  generally 
produce  a  transverse  angular  rate  of  the  rocket  as  a  whole 
(gross  motion)  subsequent  to  EOG,  the  experimental  confir¬ 
mation  of  this  phenomenon  is  currently  based  on  optical  lever 
data  acquired  without  a  priori  consideration  of  transverse 
bending  as  a  possible  error  source.  It  is  recommended  that 
a  test  program  be  conducted  to  verify  that  transverse  bending 
vibrations  are  as  significant  as  theoretically  predicted  and 
to  determine  how  well  the  current  theories  predict  their  effects. 
Such  a  test  program  would  necessarily  involve  (1)  tests  to 
determine  the  effects  of  transverse  vibration  and  (2)  tests  to 
determine  its  major  causes.  The  former  could  be  conducted  with 
small  reusable  experimental  rocket  models  which  are  propelled 
by  cold  gas,  or  spring  mechanisms,  and  which  are  caused  to 
vibrate  in  prescribed  fashions.  The  latter  could  be  conducted 
with  similar  experimental  rocket  models  and  a  launcher  for 
which  friction  and  support  stiffness  could  be  controlled. 
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APPENDIX  A 


TWO-BODY  MATHEMATICAL  MODEL  FOR  FREE  FLIGHT 

Introduction 

The  mathematical  model  developed  on  the  basis  of  the  two-body 
physical  model  described  in  Ref.  2  and  in  the  body  of  this  report  is 
given  here  for  the  convenience  of  the  reader  and  also  for  the  purpose  of 
correcting  some  errors  which  appear  in  Ref.  2.  The  derivation  of  Ref.  2 
should  be  consulted  regarding  the  nonlinear  equations  of  motion, since 
only  the  linearized  versions  are  given  here. 

Equations  of  Motion 

The  linearized  version  of  the  equations  which  govern  the  rotational 
motion  of  the  system  of  two  bodies  can  be  expressed  in  the  form  of  four 
first-order,  scalar,- ordinary  differential  equations  and  a  first-order 
matrix  equation.  The  scalar  equations  are 

SI.  =  T  /I  ,  (A- la) 

1  s  A 

$  -  £2  ,  (A- lb) 

$  =  +  ^3  (A-lc) 

and 

0  =  S^2  c$  -  s4>  ,  (A-ld) 

where  the  ^ ,  j=l,2,3,  are  the  and  2 ^-components,  respectively,  of 

the  angular  velocity  of  Body  1  (the  "aft"  body),  Tg  is  the  spin  torque; 

IA  is  the  moment  of  inertia  of  the  undeformed  system  of  two  bodies  about 

•  • 

the  x^-axis;  $  is  the  roll  angle  of  Body  1;  f  and  0  are  the  yaw  and  pitch 
angles  of  Body  1. 
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PQ II 
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The  matrix  equation  can  be  written  as  follows: 


I  x  =  B  x  +  T  , 


where 


(A-2) 


x  =  (^2  w2  o)3  e2  e3)  ; 


(A- 3) 


I  is  a  6x6  matrix  with  non-zero  elements, 


and 


Ill  =  A2+B2+a(Vr2)2> 


Z13  =  B2  +  °  r2(W’ 


I22  ~  111’ 


Z33  =  Z44  =  B2  +  °r2  ' 


Z42  "  I24  =  I31  =  I13’ 


I55  "  Z66  =  1; 


is  a  6x6  matrix  with  non-zero  elements, 

B12  =  fil[A2  +  B2  +  °<Vr2)2 


VV  ’ 


B14  =  ni^2B2  '  Bx  +  2  O  r2(Jl1+r2)], 

B15  =  yr2FT  +  V[B2~B1  +  ar2(Jl1+r2)] 

B16  ’  VV8!  *  «2<W  ' 


B21  =  "  B12’ 


B23  "  B14 ’ 


B25  *  V 


(A-4a) 

(A-4b) 

(A-4c) 

(A-4d) 

(A-4e) 

(A-4f ) 

(A-5a) 

(A- 5b) 
(A-5c) 
(A-5d) 
(A-5e) 
(A-5f ) 
(A-5g) 


r  'TJt  n-r-r '  ■■■>  ~  v 


-r—  r  ■ 
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B26 

"  B15’ 

(A-5h) 

B32 

=  ^1[b2~b1  +  ar2(s'1+r2)]  , 

(A-5i) 

B33 

=  -c, 

(A-5j ) 

B34 

=  nL[2  B2-B1  +  2or22], 

(A- 5k) 

B35 

=  pr2FT  +  ni2lB2-B1  +  ar22]  -k 

(A-5A) 

B36 

=  ^1[b2~b1  +  or2z]  ; 

(A-5m) 

T  is  the  6x1  matrix  , 

T  =  (-Ot  F  Z  -a  F  £  0  0  0  0)T.  (A-6) 

z,  yTczTc 

In  the  above  Eq.  (A-3) ,  and  are,  respectively,  the  y^-  and 

2^-components  of  the  angular  velocity  of  Body  2  relative  to  Body  1  and 

9  and  9  are  the  transverse  angles  of  rotation  of  Body  2  relative 
2  3 

to  Body  1.  Also,  in  Eqs.  (A-4)  and  (A-5),  A1  and  A are,  respectively,  the 
axial  and  transverse  centroidal  moments  of  inertia  of  Body  1,  while  B^  and 
B2  are  those  for  Body  2.  Furthermore,  with  m^  and  m2  denoting  the  masses 
of  Bodies  1  and  2,  respectively,  0  =  m^m^/M  where  M  =  m^+m2.  Additionally, 
is  the  distance  from  the  center  of  mass  of  Body  1  to  the  "hinge  point" 
between  the  bodies,  and  is  the  distance  from  the  hinge  point  to  the  center 
of  mass  of  Body  2.  Finally,  c  and  k  are  the  inter-body  damping  and  stiffness 
constants . 

In  Eq.  (A-6),  and  a ^  are  thrust  misalignment  angles;  F^  is  the 
thrust  magnitude  and  =  [m^r^  +(&^+r^+r2)m2]/M,  where  r^  is  the  distance 
from  the  aft  end  of  Body  1  to  its  center  of  mass. 
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Constant  Spin  Rate  Solution 

When  the  spin  torque  T  =  0,  fi  is  constant.  In  this  case,  if  F  is 

Si  1 

also  a  known  function  of  time,  a  closed-form  solution  to  Eqs.  (A-l)  and 
(A-2)  can  be  obtained.  Since  FT  is  fairly  constant  immediately 
after  "end  of  guidance"  (EOG) ,  it  is  assumed  to  be  constant  herein,  there¬ 
by  simplifying  the  process  of  obtaining  a  particular  solution  to  Eq.  (A-2) . 
In  the  case  of  constant  spin  rate,  B  is  constant.  Hence,  the  matrix 

A  =  I-1B  ,  (A- 7) 

is  also  constant  and  Eq.  (A-2)  can  be  written  as 

x  =  A  x  +  T  ,  (A- 8) 

where  I  =  I-1T  . 

A  particular  solution  to  Eq.  (A-8) (t  is  constant)  is 

Xp  =  ~  A  1T  ,  (A-9a) 

or 

x  =  B-1!  .  (A- 9b) 

-P  =  - 

The  complete  general  solution  to  Eq.  (A-8)  consists  of  xp  and  the  solution 
to  the  homogeneous  equation  , 

\  =  =  ?h  ’  -h(0)  =  -h0  •  (A- 10) 

By  using  the  eigenvalues  and  eigenvectors  of  A,  which  for  our  purposes 
appear  as  complex  conjugate  pairs  (three  of  them),  the  solution  x^  can  be 
expressed  as 

5h  *  I  5  !"'  5h0  '  <A-U) 

where  E  is  a  6x6  matrix  the  columns  of  which  are  alternately  the  real  and 
imaginary  parts  of  the  eigenvectors  corresponding  to  the  eigenvalues  of  A 
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with  positive  imaginary  parts.  Also,  if  n.  and  id.  denote  the  real  and 

J  J 

imaginary  parts  respectively  of  an  eigenvalue  A^  which  has  a  positive 
imaginary  part. 


2u  ^  2 

2  =  2  222  2 

2  2  £33 


(A- 12) 


where 


n  t  n.x 

e  J  cos  U).T  e  -1  gin  co  T 
J  j 

D.  .  = 

=JJ  n  t  n.x 

-e  -1  sin  cj.t  e  3  cos  u  T 

L  j  j 


(j-1,2,3) 


(A-13) 


The  solutions  for  the  first  two  elements  of  x;  i.e.,  x.,  j=l,2,  can 
be  written  out  explicitly  and  substituted  intoEqs.  (lc)  and  (Id)  to  obtain 


-7  l 

i=l 


cos  p.T  -sin  p.T  AirBi2  n  T 


[_sin  p^T  cos  pixj  |Ai2+Bi 


3  cos  q±T  -sin  q  r  A  +B 
+  7l  e  1  + 


cos<I>  -sin$  lx. 


sin  q^T  cos  q^  Ai2“Bil 


sin$  cos$  x„ 


(A- 14) 


where  $  =  T  «  t-^^  =  ^  "  t*>±,  qi  -  ^  and  if  Ejk  is  the  (j,k) 

element  of  E  and  is  the  jth  element  of  C  =  E~ 1 xQ , 


Aij  =  Ej ,2i-i  C2i-1  +  Ej ,2i  °21 


(A- 15a) 
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and 


Bij  *  Ej ,2i-l  °21  “  Ej,2i  °2i-l 


(A-15b) 


Equation  (A-14)  can  be  integrated. 


The  result  is 


"o' 

1  3  eV 

-  y  e 

P(T) 

'Ail"Bi2 

4. 

®0 

4' 

2  2^2 
1=1  n.+p. 

1  1 

W 

Jo 

~  n  T 

+  i 

g(x) 

Ail+Bi2' 

sin<l>  (cos$-l) 

Vni 

2  i=l  nj+qj 

fi2-Bil_ 

(l-cos<{>)  sin$ 

_  . 

W\ 

(A- 16) 


where 


ni^CpiT  -1)  +  pispiT 


-tnispiT 


Pi  (cp^-1)] 


P(T)  = 


(A-17) 


^  spiT-pi(cpiT-l)  ni(cpit  -1)  +  p-Sp^  | 

and  Q(t)  has  the  same  form  as  P(t),  but  in  Q(f)  the  take  the  place  of  p^. 

The  pitch  and  yaw  angles  of  the  rocket's  nose,  ©N  and  4^,  respectively, 
are  given ;in  the  linear  approximation by 


The  initial  conditions  on  the  x^ , 
obtained  by  assuming  that  the  rocket  is 
P  and  Q  which  do  not  move  transversely. 


(A- 18) 


j*l,2,...,6,  and/or  0  and  4'  are 
,  at  EOG,  supported  at  two  points 
This  constraint  and  the  assumption 


of  small  angles  lead  to  the  following  initial  conditions: 
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X10  = 

(A-19a) 

x20  = 

“[  (^2+r2^^1+rl^^NO  ’ 

(A- 19b) 

X30  = 

i««W1[%  +  Vk0i  • 

(A-19c) 

X40  = 

[  M/Uj+rp  ]  [t^o  *  ^ieno^  ’ 

(A-19d) 

X  38 

50 

[d'(lliVls»o  • 

(A-19e) 

x  = 

60 

[d/<2i+ri)]yi)0  . 

(A-19f) 

0O  '  "[(S'2+r2)/d]  X50  (A-19g) 

and 

H'q  [  (^2+r2>/d]  x60  • 
where  d  =  +  SL^  +  . 


(A-19h) 


APPENDIX  B 


ASSUMED-MODES  MATHEMATICAL  MODEL 


Introduction 

The  mathematical  model, developed  on  the  basis  of  the  "continuous" 
physical  model  and  referred  to  as  the  "assumed-mode"  model  because  of  the 
use  of  assumed  mode  shapes  to  model  the  rocket  bending, is  given  here  for 
the  sake  of  completeness.  For  additional  details,  the  reader  is  referred 
to  Refs.  2  and  3. 

Actually,  two  mathematical  models  of  a  flexible  rocket,  one  for  on- 
launcher  motion  and  a  second  for  off-launcher  motion,  are  needed.  For 
this  reason,  two  sets  of  equations  of  motion  as  well  as  the  transformation 
relating  the  two  are  given  here. 


Equations  for  On-Launcher  Motion 

When  the  rocket  is  on  the  launcher,  i.e,,  during  the  detent  and 
guidance  phases,  pinned-pinned-free  mode  shapes,  i|;  (x) ,  j=l,2,3,  are  used 
to  approximate  its  bent  shape  by  assuming  that  the  y-  and  z-displacements 
(rotating  with  the  rocket)  of  the  rocket  centerline  are 


y  =  n  +  y0(x) 

and 

z  -  S  +  zQ(x), 


(B-la) 

(B-lb) 


where 

3  3 

1  -  l  P,  (t)  (x)  and  £  *  l  q  (t)  ip  (x) 

j-1  J  J  J-l  3  3 

are  the  parts  of  y  and  z,  respectively,  due  to  bending  of  the  rocket  and 


yQ(x)  and  zQ(x)  are  functions  which  define  its  original  non-straight 
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shape  if  necessary.  The  p_.  and  q^  are  generalized  coordinates  which  are 
six  of  the  sixteen  (16)  state  variables  which  are  used  to  model  the 
rocket's  motion. 

The  equations  governing  axial  rotation  are 


ft.  -  T  /I 
1  s  A 


(B-2a) 


* 


(B-2b) 


where  Tg  is  the  spin  torque  and  IA  is  the  axial  moments  of  inertia  of  the 
unbent  rocket.  Here,  ft^  is  the  x-component  of  the  angular  velocity  of  the 
rocket  and  $  is  the  rocket  roll  angle. 

Simple  equations  are  also  used  to  model  the  translational  dynamics  of 
the  point  P  (on  the  centerline  of  the  rocket  at  its  aft  end) .  These  are 


*P  =  UI 


(B-3a) 


Up  =  Ft/M  -  g  sin  0L, 


(B-3b) 


where  Fp  is  the  thrust  magnitude  and  0^  is  the  fixed  launcher  elevation. 

The  equations  for  the  p^  and  q^.  are  obtained  as  explained  in  Refs.  1 
and  2.  The  results  are 


M  (p  “2ft  q-ft  q-ft2p)  +  (K-FJ)p  -  (g  sin  $)m  +  F  c  -fib-  fi2a 

X—  X—  X—  ~  i—  —  “  1“  X—  X"* 


+  Dp  =  0 


(B-4a) 


M(q+2fi  p  +  fi  p  -  fi2q)  +  (K-F_  Dq 
»  -  x£  x-  x-  =  T  =  - 

-  (g  cos  <I))m  +  FTd+ftxa-fi2b+Dq=0 


(B-4b) 


T  T 

where  g  *  (p^  p^  P^)  and  q  -  (q^  q^  •  Here  also,  M  is  (currently)  a 
3x3  matrix  with  elements 
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*  /  a \p  \b  dx 

J  n  1  J 


(B-5a) 


where  a  is  the  mass  per  unit  length  of  the  rocket;  K  is  a  3x3  matrix  with 
elements 


K .  .  =  /  i|»,(E  I  dx  , 

13  ~  l  3 


(B-5b) 


where  (  )'  =  d(  )/dx  and  El  is  the  bending  stiffness  of  the  rocket  at 
station  x;  T  is  a  3x3  matrix  with  elements 


XJ 

y . .  =  /’/>!<('!  f(x)  dx  , 

n  1  -J 


(B-5c) 


where  f(x)  is  the  internal  force  per  unit  length;  m  is  a  3x1  matrix  with 


elements 


m  =  /  a  ip.  dx 

1  0  1  • 


(B-5d) 


Furthermore,  the  matrices  a,  b,  c  and  d  are  3x1  matrices  with  elements, 


h  m  I  ^  <|j  y  dx  , 
0 


(B-6a) 


bi'  Js  h  zo • 


(B-6b) 


ci  ■  J  5  *1  yo 


(B-6c) 


di  3  /  5  ip'  z’  dx  , 

0 


(B-6d) 


respectively.  Also,  g  is  the  gravitational  acceleration  magnitude. 

The  angular  rotation  of  the  nose  of  the  rocket  is  defined  by  the 
angles  0^  and  'F  as  in  the  case  of  the  two-body  model.  However,  for  the 
continuous  model  these  angles  are  determined  by  using 
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0  =0+0  cos  $  -  0  sin  $ 

N  L  2N  3N 


(B-7a) 


\  =  92N  Sln  *  +  03N  C°S  *  ' 


(B-7b) 


where 


e2N  =  qi  K(L)  ~  Z0(L) 


(B-8a) 


93N  '  I  »i  *ia)  +  • 


(B-8b) 


Equations  for  Off-Launcher  Motion 

The  motion  of  the  rocket  after  end  of  guidance  is  modeled  by  using 
mode  shapes  for  a  free- free  uniform  beam  and  using  instantaneous  principal 
axes  of  the  deformed  rocket  in  writing  the  equations  for  rotational  motion. 
The  use  of  instantaneous  principal  axes  simplifies  the  equation  of  rotational 


motion  so  that  it  is  simply 


Ic  0)  +  w  Ic  0)  =  T  , 


(B-9) 


where  I  is  the  principal  centroidal  inertia  matrix,  w  =  (w  to  oj  )T  is  the 
-C  i  2  3 

matrix  of  u^-,  u2~  and  u^-components  of  the  angular  velocity  of  the  principal 


(Cu-jU^)  system, 

0  -u>3  co2 

<5  =  0  -w 

-  3  1 

-<*>2  0 


(B-10) 


T  =  {Ts  Ftxc  lu'(O.t)  +  u3(xM,t)/xc  -  a  ] 


-  Ft  xc  [u^(0,t)  +  u2(xM,t)/xc  +  aj}1  . 


(B-ll) 
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In  Eq.  (B-ll) ,  x  is  the  distance  from  point  P  to  the  center  of  mass  of 
the  rocket,  x^  is  the  distance  from  P  to  the  head  end  of  the  rocket 
motor,  u^(0,t)  andu^(0,t)  are  the  slope  of  the  bent  (and  bending) 
centerline  of  the  rocket  at  P  and  a  and  a  are  constant  mechanical 

y  z 

thrust  misalignment  angles. 

The  translational  motion  of  the  rocket  is  governed  by  the  equation, 


v=-iv+  F  /M  +  g  ,  (B-12) 

where  V  is  the  3x1  matrix  of  u^-  ,  u^-  and  u^-components  of  the  velocity 
of  the  center  of  mass  C,  g  is  a  3x1  matrix  of  similar  components  of  the 
acceleration  of  gravity  and 


Ft  =  F  [  1  u'(0,t)+a  u' (0, t)  -  a  ]T.  (B-13) 

-l  i  j  z  l  y 

In  addition  to  Eqs.  (B-9)  and  (B-12),  the  following  kinematic 
equations  are  needed-  to  complete  the  gross  motion  equations: 

Rotation 


+  (uj2  sin  <J>  +  ojj  cos  <p)  tan  0 


9  =  w2  cos  4>  —  oo^  sin  <p 

=  (oj^  sin  4>  +  <a  cos  <J>)/ cos  0 


Translation 


cip 

-sip 

Sip 

cip 

0 

0 

c0  0  s9 

0  10 
-s9  0  c0 


As  in  the  case  of  on-launcher  motion. 


(B-14a) 

(B-14b) 

(B-14c) 


m' 

o 

o 

_ 1 

u 

0  c^>  -s<j> 

V 

0  s<{>  c<jj 

w 

it  is  assumed  that  the 


(B-15) 


rocket 


may  be  non-straight  before  any  bending  occurs.  Hence,  the  functions 
u2(x,t)  and  u^(x,t)  are  expressed  as  follows: 
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and 


u2(x,t) 


3 

=  u2Q(x)  +  l  P-  (0<l>.  (x) 


3 

u3(x,t)  =  u3Q(x)  +  l  Qj  (x)  , 


where  u2q(x)  and  u3q(x)  are  the  counterparts  of  y^Cx)  and 
Pj  and  Qj  are  new  generalized  coordinates,  i.e.,  they  are 
such  that  P^  =  p^,  etc. 

The  ways  in  which  the  P^  and  change  with  time  are 
solving  the  equations , 

=  f  *  2W!  9  “  “j.  9  ~  <“ i  P)  +  (Kf  -  ft  T)P  +  DP 

+  (oi3  +  ~  kf  ~  5f  +  ft  £f 

-  F_  [u'  (0,  t)  +  a  1  <f>(0)  +  (F_/M)_  m  =  0 

1  L  Z  “  1  Z  —I 


and 


Mf  (Q  +  2w.jP  +  w1  P  -  w12Q)  +  (Kf  -  FT£f)g  +  DQ 
+  (“3  -  “2>-f  +  wi  kf  “  Vkf  +  ft  -f 

+  FT[u^(0,t)  +  o  ]  +  (Ft/M)3  mf  =  0, 


where  a  subscript  f  is  used  to  denote  "f ree-f light. " 

The  matrices  appearing  in  Eqs.  (B-17)  are  defined  as 
P  and  Q  are  3x1 ,  i.e., 

.  T 


and 


P  =  (Pj_  ?2  P3) 


9  =  (Qx  Q2  Q3) 


(B-16a) 

(B-16b) 

z^(x)  and  the 
not  necessarily 

determined  by 


(B-17a) 

(B-17b) 


follows : 

(B-18a) 


(B-18b) 
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M.  is  3x3,  with  elements. 


m .  ,  =  /  a  .  <fc .  dx  , 

1-4  J  -t  -5 


‘f  0 


i  J 


is  3x3,  with  elements, 

L 

K..  =  /  <+>  [El  dx, 

1Jf  0  1  J 


£  is  3x3,  with  elements. 


Yii  =  ^  dx  , 


Jf  0  1  7 


a^,  bf,  c^,  d_,  and  are  3x1,  with  elements. 


a  =  /  5  <j»  u  dx  , 

f  0 


blf  '  J  5  *1  U30  dX 


1  '  J  3  *1  “2 


0 


20 


dx 


and 


di£  ‘  J  5  *i  So dx 


i  =  f  a  <t>  (x-x  )  dx 
Xf  0  1  L 


m  =  f  a  <p  dx  , 
f  0 


respectively . 


(B-18c) 

(B-18d) 

(B-I8e) 

(B-19a) 

(B-19b) 

(B-19c) 

(B-19d) 

(B-19e) 

(B-19f ) 
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Also  in  Eqs.(B-17), 

(F_/M),  =  (F/M)[u  (0,t)  +  a  ]  (B-20a) 

T  Z  T  i.  Z 

and 

(Ft/M)3  =  (FT/M)[u3(0,t)  -  a  ]  .  (B-20b) 

Transformation  Required  at  End  of  Guidance 

Since  two  sets  of  equations  of  motion  are  used,  one  for  the  guidance 

phase  and  one  for  the  free-f light  phase,  it  is  necessary  that  the  solutions 

to  these  equations  be  matched  at  the  EOG.  Hence,  we  require  that  the 

positions  and  velocities  of  points  P  and  Q  and  of  the  center  of  mass  be 

continuous  at  t^.^,  the  time  of  EOG.  Since  three  mode  shapes  are  used 
EOG 

during  each  phase,  we  can  also  require  that  the  angular  orientation  and 

angular  velocity  of  the  nose  of  the  rocket  be  continuous  at  t„^_. 

To  enforce  the  desired  continuity  of  state,  we  must  first  determine 

the  orientation  and  angular  velocity  of  the  centroidal  principal  axes  of 

the  rocket,  because  the  bending  during  free-f light  is  referenced  to  those 

axes.  The  centroidal  principal  axes  are  oriented  with  respect  to  the  Pxyz 

system  by  using  the  angles  an(*  ^3  which  are  defined  as  follows: 

L  _ 

P2  *  -[  /  a(z0+S-zc) (x-xc)dx]/IT  (B-21a) 

and 

L  _ 

y3  =  [  /  5(yQ+n-yc) (x-xc)dx]/IT  ,  (B-21b) 

L  L  L 

where  x  =  T  /  axdx]/M,  y  -  [  /  a(y  +n)dx]/M  and  z  =  C /  a(zn+Odx]/M. 

L  0  0  U  L0 

The  orientation  of  the  Cu^u^u^  system  in  space  is  defined  by  the 

angles  ip,  Q  and  <J)(3-2-l  sequence).  At  EOG,  these  angles  can  be  expressed 

in  terms  of  (the  roll  angle  of  the  Pxyz  system)  an(*  by  equating 
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elements  of  equivalent  direction  cosine  matrices.  The  required 
expressions  are 


4>  =  tan-1  ((U3c4>  +  u2s$)  /  [cGM-Cu^  -V^cfc)  sGJ  }  , 


9  =  sin-1[s6^  +(u2c$  -  y3s$)c0^] 


(B-^2-0 

(B-22b) 


<p  =  tan-1  {(-y  s0  +  s4>c0  )/(-y  s0  +  c4>c9  ) } 
0  Li  L  Z  Li  Li 


(B-22c) 


Since  the  angular  velocity  of  the  cuj_U2U3  system  is  that  of  the  Pxyz 
system  plus  the  angular  velocity  of  the  cu^u^u^  system  relative  to  the 
Pxyz  system, 

<*>  ~  t  U2  ~  U3^x  U3  +  U2ax]T  •  (B-23) 

The  position  and  velocity  of  the  center  of  mass  C  are  obtained  by 


using  the  definitions  of  x  ,  y  and  z  given  above.  The  space-fixed 

L>  C  C 

coordinates  X  ,Y  and  Z  of  C  while  the  rocket  is  guided  are  found  from 

L*  L*  L» 

the  equation, 

fxT  f  c0_  0  s0l  fl  0  Ol  Ix.+Xl 


c0  0  s0,  1  0  O'  Tx^+X  I 

L  L  C  p 

0  10  0  c$>  -s<J>  Y_ 


s0T  0  c0  0  s<f  c$  Z 
L  L  C 


(B-24) 


Furthermore  the  u^-components  of  the  velocity  of  C  are,  during  guidance, 


u  =  *P 


v  =  -y3u  +  yc  -  Q^zc 


(B-25a) 

(B-25b) 


w  =  y.u  +  z  +  n  y_. 

Z  L  X  0 


(B-25c) 


Because  only  a  limited  number  of  mode  shapes  are  used,  the  position 
and  velocity  of  each  particle  of  the  rocket  cannot  be  forced  to  be 
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continuous  at  EOG.  However,  if  three  mode  shapes  are  used,  the  position 
and  velocities  of  three  points  on  the  rocket's  centerline  can  be  made 
continuous  at  EOG.  Alternatively,  the  positions  and  velocities  of  two 
points  on  the  rocket's  centerline  and  the  slope  and  time  rate  of  change 
of  slope  of  a  third  point  can  be  made  continuous  at  EOG.  The  latter 
approach  is  taken  because  it  is  desired  that  the  pitch  and  yaw  angles  of 
the  nose  of  the  rocket  and  their  time  rates  of  change  be  continuous  at 
EOG.  Also,  to  model  the  effects  of  the  constraints  during  guidance,  the 
positions  of  points  P  and  Q  are  forced  to  be  continuous  at  EOG. 

If  we  let  Xq  be  the  position  of  the  forward  support  and  def:  the 

matrices 


■*,(0) 

4>2  (0) 

4>3(0) 

$  = 

*2(V 

<03(Xq) 

<^(L) 

^•(L) 

— 

~i 

b  * 

-z 


-y= 


(B-26a) 


"M3Xc 

(VXC)U3"yc 

p  +  n' (L,t)  +  y' (l) 

->  t=t 


(B-26b) 


EOG 


Vc 

-U2  (VV  'ZC 

l-U,  +  (L,t)  +  z' (L) 

L  U 


(B-26c) 


EOG 


1 

(B-26d) 
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and 


b  = 

-z 


->'2(VtC)  “c 

-P2  +  C'(L,t) 


t=t 


EOG 


then. 


5<cEOG>  ■  r‘  iy  ' 


9<w>  -  r‘  6Z . 


EOG 


p(t  )  =  $-1  b 
-  EOG  =  -y 


and 


Q(t  )  =  b  . 
-  EOG  =  -z 


(B-26e) 


(B-27a) 

i 

(B-27b) 

(B-27c) 

(B-27d) 


The  above  "initial"  conditions  (those  at  EOG)  on  X  ,  Y  ,  Z  ,  u,  v,  w, 

•  c  c  c 

ip,  9,  <j>,  u),  P  and  Q  can  be  used  to  continue  the  solution  of  the  equations 
of  motion  during  the  free-f light  phase. 


APPENDIX  C 


PRINCIPAL  AXES  ROTATIONS 
DUE  TO  TRANSVERSE 
VIBRATION 

Introduction 

Due  to  slenderness  of  most  free-f light  rockets,  there  is  generally 
little  coupling  of  their  transverse  vibrational  and  rotational  motions 
subsequent  to  end-of-guidance  (EOG)  if  aerodynamic  forces  are  not  signif¬ 
icant.  Hence,  as  explained  in  the  body  of  this  report,  the  free-f light 
motion  of  a  slender  free-flight  rocket  which  is  caused  by  bending  vibrations 
prior  to  EOG  is  essentially  composed  of  a  uniform,  transverse  rotational 
motion  of  the  rocket  about  its  centroidal  principal  axes,  a  vibrational 
motion  with  respect  to  those  axes,  and  a  translation  of  the  rocket's 
center  of  mass.  In  this  appendix,  we  present  an  approximate  method  for 
determining  the  transverse  angular  velocity  of  the  rocket's  principal  axes 
at  EOG  and,  hence,  the  transverse  rate  due  to  transverse  vibration.  The 
two-body  model  described  in  Ref.  2  and  herein  (see  Section  2  and  Appendix 
A)  is  used. 

Principal  Axes  of  the  Two-Body  Model 

The  orientation  of  the  centroidal  principal  (see  Fig.  C-l)  reference 
frame  Cx'y'z'  relative  to  the  Cxyz  reference  frame  (which  is  parallel  to 
the  Body  1-fixed  C^x^y^z^  reference  frame)  can  be  defined  in  terms  of  the 
inertia  characteristics  of  the  bodies  and  the  angles  @2  and  0^. 
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6 


Fig.  C-l.  Principal  axes  of  the  two-body  model. 

The  inertia  matrix  of  the  system  of  two  bodies  using  the  Cxyz  frame 
as  the  reference  is, through  first  order  in  $2  and  0 


where 


lA  =  A1  +  B1  * 


the  sum  of  the  axial  moments  of  inertia  of  Bodies  1  and  2,  and 


J  =  B2  -  Bi  +  o  r2^1+r2)  • 


(C-l) 


(C-2) 


(C— 3) 


In  Eq.  (C-3) ,  B2  is  the  centroidal  transverse  moment  of  inertia  of  Body  2 
cr  =  m1m2/(m1-Hn2)  where  m.  is  the  mass  of  Body  j,  r2  =  |r2|  and  A,  -  |£j. 
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Furthermore, 


XT  =  A2  +  B2  +  a<£1+r2)2’  (C-4) 

where  is  the  centroidal  transverse  moment  of  inertia  of  Body  1. 

Now,  we  let  I'  denote  the  principal  centroidal  inertia  matrix  of  the 
system  and  and  denote  small  angles  which  define  the  orientation  of 
the  centroidal  principal  frame  Cx'y'z'  with  respect  to  the  Cxyz  frame. 
Then,  by  using  the  fact  that 


1  ~U3  U2‘ 

1  U3  V 

I  _ 

3  — 

u3  1  0 

IIH 

-U3  1  0 

_-U2  0  1 

— 1 

O 

CM 

_ 1 

we  get, through  first  order  in  and  y^, 

U2  “  tJ/(VIA)]62  (C-5a) 

and 

W3  =  tJ/(W103  *  (C“5b) 

Angular  Velocity  of  the  Principal  Frame 

The  angular  velocity  of  the  principal  frame  Cx'y'z',  as  resolved  in 
terms  of  components  in  that  frame,  is 

U) '  = 

where  an<*  ^3  are  the  x-,  y-  and  z-components,  respectively,  of  the 

angular  velocity  of  the  Cxyz  reference  frame.  By  neglecting  terms  which 

•  • 

are  second  order  in  y3*  y2  anc*  ^3*  we  8et  the  following  expressions 
for  the  components  of  u>' : 


1  U3  "U2 

^7 

1 _ 

0 

-y3  1  0 

+ 

• 

U2 

y2  0  1 

_^3 
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-v 

(C-7a) 

“2 

‘  S2  +  “2 

(C-7b) 

W3 

=  +  u2^1  +  u3  . 

(C-7c) 

Equations  (7)  may  be  used  at  any  time,  provided  the  assumptions  of 

small  angles  and  angular  rates  made  above  are  valid.  We  are,  however,  most  J 

interested  in  the  angular  velocity  of  the  principal  axes  reference 

frame  at  the  time  of  EOG,  since  if  there  is  no  thrust  misalignment  and 

negligible  aerodynamic  torque,  the  angular  rates  in  pitch  and  yaw  present 

at  EOG  will  persist  during  free-f light. 

If  the  constraints  adopted  are  such  that  points P  and  Q  (see  Fig.  Al) 

can  only  translate  along  a  fixed  straight  line  during  the  guidance  phase, 

then  at  the  EOG  when  t=t„_„  (and  letting  a  subscript  0  denote  the  value  of 

EOG 

a  variable  at  that  time) 


(C-8a) 
(C-8b) 
(C-8c) 
(C-8d) 
(C-8e) 
(C-8f ) 
(C-8g) 

(C-8h) 
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In  Eqs .  (8) , 


and 


where 


and 


ki  ■  d/di 


k2  =  d/d2> 


d  =  dl  +  d2  ’ 
d!  "  41  +  ri 


d2  '  *2  +  r  • 


(C-9a) 

(C-9b) 

(C-lOa) 

(C-lOb) 

(C-lOc) 


From  Eqs.  (C-5)  and  Eqs.  (C-8)  we  find  that  at  EOG,  the  values  of 


y2  and  y^  are 


and 


d20  ‘  S  kl  ®SO 


30 


5  kl  ^NO* 


respectively,  where 


5  =  J/(It-Ia)  , 


and  the  values  of  y2  and  y^  are 


U20  ~  6  kl'ni  \o  + 


(C-lla) 


(C-llb) 


(C-12) 


(C-13a) 


and 


1*30  ■  «  V®!  0„o  +  »„>  , 


(C-13b) 

respectively.  By  using  Eqs.  (C-ll)  and  (C-13)  in  Eqs.  (  C— 7  ),  we  get  the 
following  expressions  for  the  values  of  and  at  EOG: 

(C-14a) 
(C-14b) 


“20  '  <6kl  -  kl/k2)9NO 


“30  ‘  (Sk!  -  VW 


Note:  Here,  &2  such  that  £2  *  ^2  d2*  i*e*i  it  is  the  X2~component 
of  the  vector  from  C2  to  Q.  Hence,  it  may  be  negative. 
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Also 


The  transverse  angular  rate  due  to  transverse  vibrations  is  therefore 


“t  ‘  ’so'  <C-15> 

the  aim  change  due  to  transverse  vibration  is 

0„  =  <5k\/e2  +  f2  /r  ,n 

t  r*  NO  NO  (C-16) 


APPENDIX  D 


TWO-BODY  MATHEMATICAL  MODEL 
WITH  FLEXIBLE  SUPPORTS 


Introduction 

Since  the  structure  which  supports  a  rocket  until  end  of  guidance 
(EOG)  is  in  reality  not  rigid,  the  issue  of  whether  the  flexibility  of 
this  structure  contributes  to  the  bending  of  the  supported  rocket  needs 
to  be  addressed.  Furthermore,  the  effects  of  frictional  forces,  which 
necessarily  act  on  the  rocket  at  its  points  of  support,  must  be  considered 
as  a  possible  mechanism  for  causing  transverse  vibration  of  the  rocket. 

To  this  end,  the  two-body  model  presented  in  Appendix  A  was  extended  to 
include  the  effects  of  support  flexibility  and  friction. 

Extended  Model 

The  physical  model  is  depicted  in  Fig.  D-l.  It  is  identical  with  that 
used  in  deriving  the  equations  given  in  Appendix  A,  except  the  rocket  is 
supported  in  such  a  manner  that  forces  are  exerted  on  it  through  small 
"shoes."  The  shoes  in  the  present  model  are  not  constrained  to  move  on 
rails,  but  are  assumed  to  slide  on  the  interior  of  a  "smooth"  tube. 

Although  the  tube  is  smooth,  it  is  not  perfectly  so;  hence,  frictional 
forces  are  assumed  present  at  the  rocket/tube  interface. 

Equations  of  Motion 

Because  the  points  P  and  Q  can  move  transversely  as  well  as  longitudi¬ 
nally,  a  general  equation  for  motion  of  the  mass  centers  of  the  two-body 


Fig.  D.l.  Two-body  model  in  flexible  tube, 
system  is  required.  The  center  of  mass  C  is  located  by  the  vector, 

Rc  =  Rp  =  yai+r2)  +  £1,  (D-D 

where  y  =  n^/Giyfn^)  and  m  is  the  mass  of  Body  j.  Here  (as  in  Appendix  A), 
the  aft  body  is  Boav  1.  The  acceleration  of  C  is 


Mr,  1 

2du) 

h  =  h  +  u 

la 

--  x  (V£i)  +  5X^(-1+I1)  ] r 

u  —  X  £2  + 

(D-2) 


where  J2  is  the  angular  velocity  of  Body  1,  w  is  the  angular  velocity  of 
Body  2  and  ^d(  )/dt  denotes  the  derivative  with  respect  to  time  of  (  ) , 
treating  the  unit  vectors  associated  with  Cj y j z ^  as  non-rotating. 


i 
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If  we  let 


A  = 


1 

-03  1 


02  0 


0-  -9, 


and 


0 

1 


(D-3) 


l 

V 

-0 

_C0L 

0 

-S0l" 

£  = 

-fc^H-Gs^ 

C  $ 

s$ 

0 

1 

0 

-s$ 

c  $ 

_s0L 

0 

C0L. 

(D-4) 


denote  matrices  which  define  transformations  from  the  Body-1  basis  (i  ,j  ,k  ) 

/V  /V  /V  —1  1  —1 

to  the  Body-2  basis  (A2’i2’-2^  an<^  ^rom  t^e  fixed  basis  (I,J,K)  to  the 


A  /\ 


launcher  basis  (i^,j^,k^),  respectively,  then  Eq.  (D-2)  can  be  replaced  by 
the  matrix  equation  (launcher  basis). 


?p=  ~C^  ytftftf^  +  r^  +  4^£2-2/l  =  -  ~^l+-T  ?2-^ 

-at  i2  “2/i]  +  -p/M  +  IQ/M  +  ?t/m  +  §  +  2T  £x  0  £T  [%!»  (D_5) 

where  ft  =  (ft;L  ft2  ft3>T,  0  0)T,  u  =  w  +  Aft,  u>  =  (0  92 

r2  =  (r2  0  0)T,  Fp  =  (F  Fp  Fp  )T,  F  =  (F  F  F  )T, 

*L  yL  ZL  \  \  \ 

-T  =  (FT  azFT  _ayFT^T’  M  =  mi+m2  and  i  =  8  sin  ®l  0  8  COS0L)T- 

Here,  and  are  thrust  misalignment  angles.  Also, 


Rp-(xp  5p  6p  )\ 


P  P  P 

yL  ZL 


(D-6) 


where  6  and  5  are  the  displacements  of  P  in  the  y  -  and  z  •  directions, 

yL  ZL 

respectively. 
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Equation  (D-5)  can  be  used  to  find  5p  and  5^  ,  if  Fp,  Fg,  F^,  ^/x  anc* 

yL  \ 

are  expressed  in  terms  of  the  state  variables,  >  6p  ,6p  , 

•  •  •  •  .  P  yL  ZL 

f^,  fi2’  n3>  w2  =  92’  ^3  =  63’  XP‘  5P  and  6P  and  time'  We  Sha11 

yL  \ 


consider  Fp  and  Fq  later. 

The  equations  for  the  rotational  motion  of  the  system  of  two  bodies 
about  their  center  of  mass  are  the  same  as  in  the  "free"  case  considered  in 
Appendix  A,  except  an  additional  torque  is  now  applied  to  the  two-body 
system.  It  is  due  to  the  forces  on  the  system  at  P  and  Q  which  produce 
transverse  and  axial  (x^-axis)  torque  components  and  can  be  expressed 
as 


-Tsysadd-  <-°l  '  Il>  *  -P  +<£2+-V  *  +  Ip  +  Iq  '  <D-7) 

where  Tp  and  Tg  are  torques  about  P  and  Q,  respectively.  They  are  assumed 
to  be  axial  and  due  to  friction. 

Additional  terms  must  also  be  added  to  the  equations  which  govern  the 
rotation  of  Body  2  relative  to  Body  1.  These  can  be  found  by  considering 
the  equation  for  rotational  motion  of  Body  2  about  the  accelerating  point  R 
We  have 


»2/R  '  Ir  *  "2  H2  *  iR  • 


CD-8) 


where  is  the  time  rate  of  change  of  the  angular  momentum  of  Body  2  due  to 

rotation  about  R,  T^  is  the  torque  about  R,  and  a^  is  the  absolute  accelera¬ 
tion  of  R.  The  acceleration  a^^  can  be  expressed  in  the  form, 


2r  =  h  +  +  h  • 


(D-9) 
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Also, 


-R  =  -2/1  +  -Q  +  (H2+^2)  X  Iq  +  m2-2  X  I  * 


(D-10) 


where  ^2/1  c^e  tor<?ue  on  Body  2  due  to  the  torsional  springs  and 
dampers  used  to  represent  the  flexibility  of  the  rocket  being  modeled. 
Now,  the  acceleration  Rp  can  be  expressed  as  [see  Eq.  (D-l) ] , 


Rp  =  Rc  -  fi  -  Oj+r^  ,  (D-ll) 

where  y  =  m^/M.  Obviously, 

Sr-Sc-  <Vi2>-  (IK12) 

However,  we  also  have,  from  Newton's  second  law, 

5C  =  §  +  (-T  +  -P  +  -q)/m’  (d_13) 

and  by  using  this  result  along  with  Eq.  (D-10)  in  Eq.  (D-8) ,  we  get 

§2/r  ■  I2/1  +  I,  -  >*  £2  *<?t  +  Ep>  +  l(1-u)E2  +  -2^  1  !q  •  <D-14) 

The  additional  torque  on  Body  2  due  to  the  flexible  supports  and  friction 
is  therefore 


T  =  T  -  yr  x(F  +F  )  +  [(l-y)r  +i.  ]  x  F  .  (D-15) 

add  ~  _  ~ 

Definition  of  the  Forces  and  Torques  at  P  and  Q 

The  forces  at  P  and  Q  are  assumed  to  be  caused  by  deformation  of  the 

supporting  structure,  e.g.,  tube,  and  coulomb  friction.  The  possibility 

that  the  supporting  structure  will  not  always  be  in  contact  with  the 

rocket  at  either,  or  both,  points  P  and  Q  is  allowed  for  by  simply  enclosing 

the  rocket  in  a  deadspace  "cylinder"  of  thickness  e.  For  the  purpose  of 

illustrating  how  the  forces  at  P  and  Q  are  defined,  figures  which  imply  the 
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rocket  is  within  a  tube  are  used.  However,  the  model  of  the  forces  is 
also  representative  of  other  types  of  launchers. 

In  Fig.  D-2,  the  displacement  of  point  P  from  the  launch  axis 

A  A 

(x^-axis)  is  represented  by  6p  and  6p  in  the  (y^-axis)  and  k^z^-axis) 

ZL 

directions,  respectively.  Let 


U  ■  5p  k +  {p  -l  ■ 
H  ZL 


(D-16) 


We  say  that  if 


A 

'iP  •  il >e’ 


(D-17) 


where  j  is  a  unit  vector  which  rotates  with  the  rocket,  then  the  component 
of  the  elastic  restoring  force  in  the  y-direction  is 


N  =  -  K  5  •  j  , 

Py  P  -P  - 


(D-18a) 


where  Kp  is  a  constant  in  this  study,  but  could  be  considered  a  function  of 
the  distance  traveled  along  the  tube.  Similarly,  the  z-component  of  the 
restoring  force  is 


N_  *  -  K_  6  •  k 

P  P  -P 

z 


(D-18b) 


if  j  6p  •  k |  >  e,  where  k  also  rotates  with  the  rocket. 


Rear  View 


Fig.  D.2.  Displacement  of  point  P. 
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Referring  to  Fig.  D-3,  if  $  (the  spin  rate)  is  nonzero,  frictional 

forces  which  oppose  such  spin  are  present  i^  the  rocket  is  in  contact  with 

•  • 

the  supporting  structure.  If  Up  t  0  ,  <J>  >  0  and  Up  *  0  for  <J>  *  0,  we 

have  the  transverse  frictional  force  components, 

fP  =  "  NP  (D-19b) 

y  z 

and 

fP  "  yp  NP  ‘  (D-19c) 


Fig.  D.3.  Normal  and  frictional  forces. 


In  addition  to  the  restoring  and  coulomb  friction  forces,  a  viscous 
damping  force  is  used  to  approximately  model  structural  damping.  The 
transverse  components  of  this  force  are 

dp  «  -SpKp  d(6p  •  j)/dt  0>-20a) 

y 

and 

dp  “*?pKp  d(6p  •  k)/dt  , 
z 

where  Spis  a  constant,  and  the  axial  component  is  zero. 


(D-20b) 
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If  the  rocket  is  translating  in  the  x^  direction,  then  an  additional 
force  due  to  friction  acts  at  P  in  the  negative  ^-direction.  The  x- 
component  of  this  force  is 

Fp  =  Up[  |N  I  +  |Np  |]  ,  (^  >  0).  (D-21) 

x  y  y 

The  forces  acting  at  Q  have  the  same  forms  as  those  acting  at  P;  hence, 
if  the  subscript  P  is  changed  to  Q  in  Eqs.  (D-18)  through  (D-21)  the  forces 
acting  at  Q  are  obtained.  However,  the  displacement  of  the  point  Q 


+  4'(r1+«-1+r2+«.2)  +  [03c^2s$]  (r2+«-2) 


(D-24a) 


5  -5  -  0(r1+H1+T2+«.2)  +  [63s$  -  e2d>]  (r2+^2)  (D-24b) 

ZL  ZL 

The  torques  on  Bodies  1  and  2  due  to  all  additional  (those  not 
present  in  free  flight)  sources  can  be  found  by  using  Eqs.  (D-7)  and  (D-14). 
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The  total  torque  on  the  system  about  its  center  of  mass,  obtained  assuming 
that  0,  and  0^  are  small  angles,  can  be  expressed  as 


T“Ti+Ti+Tk 
-C  A  -1  1  -1  2  -1  ’ 


(D-25) 


where 


h  *  (TP  +  TQ  +  Ts  "'Wp  +  »3?P  1  +  ][0J  ^Vn  I}  • 


2  2  2  Q  3  0  J 

y 


T1  =  {v92r2FP  +  lrl+u(5'l+r2)]Fp  ~t  d-y)r2  +  A2]02F 

- ^  Z  TC 

-[(1-U)r2  +l2lPQ  -  a  xc  PT  +  ^  rs  Sp 


(D-26a) 


+  y„  r  N 

Q  s  Q„ 


(D-26b) 


T2  *  {w93r2FP  -tr  -Hid  +r  )]Fp  -  [(l-y)r  +  «,  ]  9  F 

_ ?c  ■L  *■  y  z 

+  [(l-U)r.  +  f,]  F  -ax  F  -  y  r  N 
2  2  Q  z  C  T  P  s  P 

y  y 

-VsV  ■ 


(D-26c) 


where  Tp  and  are  axial  torque  components  due  to  friction  at  P  and  Q, 
respectively,  Tg  is  the  spin  torque,  assumed  to  be  applied  by  an  eroding 
spin  turbine, 


Np  ,  *  -  0  , 

?  ’  Y 

y  Np  -ypNp  ,  $  >  0 

y  z 


(D-27a) 


Np  ,  $  -  0  , 

a  Z 

’z  Np  +UpNp  ,  $  >0. 

z  y 


(D-27b) 


Force  components  with  Q  subscripts  are  of  the  same  form  as  those  with 
subscripts.  The  underlined  terms  inEqs.  (D-26)  are  zero  if  Xp  -  0.  They 


represent  the  "couple"  parts  of  torque  components  due  to  frictional  forces 

acting  axially  at  P  and  Q,  whereas,  the  terms  involving  Fp  and  F  represent 

x  ^x 

the  remainders  of  the  components.  Here  rg  is  the  distance  from  the  rocket 
centerline  to  the  end  of  a  "shoe." 

In  a  similar  manner  the  transverse  torque  on  Body  2  can  be  expressed  as 


T  =*  T  i  +  T  k 
-2  3  -1  4-1  ’ 

trans 


(I>-28) 


where 


T.  -  {ur.[0,Fp  +  F  1  -  [(l-u)r+L][8Fn  +F„  ] 


2  2  P  P  1 

_ X  2 

-c52  -k02  +  “qVq,1 


2  2Jl  2  Q  Q  J 
_ 2x  xz 


(D-29a) 


T4  -  {^2[03Fp  -F  ]  -  [(l-u)r  +£  ]  [0  F  -F  ] 

_ x  y  _ > 

-ce3  -  k03  Uqrs 


(D-29b) 


where  again  the  underlined  terms  are  zero  if  Xpm0.  Here,  as  in  Appendix  A, 
k  is  the  stiffness  of  the  torsional  springs  which  connect  the  two  bodies  and 
c  is  the  viscous  damping  coefficient. 

Final  Form  of  the  Mathematical  Model 


The  complete  mathematical  model  used  to  obtain  the  results  presented 
in  this  report  regarding  the  effects  of  support  flexibility  and  friction 
as  follows: 


| 
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Equations  of  Motion 


♦  =  n1. 

(D-30a) 

• 

0  3  ^2  cos<I>  -  ^3  sin 

(D-30b) 

3  ^2  sin<$  +  Q3  cos 

(D-30c) 

XP  =  UP’ 

(D-30d) 

5p  “  Cp  . 

(D-30e) 

.  yL  yL 

5P  -S  • 

(D-30f ) 

ZL  ZL 

3  T  /I  [T  given  by  Eq.  (D-26a)], 

(D-30g) 

x 3  r1  b  x  +  r1  t  , 

(D-30h) 

where 

I  3  inertia  matrix  defined  in  Appendix  A, 

B  3  coefficient  matrix  defined  in  Appendix  A, 

X  3  («2  a3  w2  u>3  92  63)T,  (D-30i) 

and 

T  3  (Tx  T2  T3  T4  0  0)T  T,  j«l,2,3,4 

defined  by  Eqs.  (D-26b) ,  (D-26c) ,  (D-29a)  and  (D-29b) ,  respectively. 


(D-30k) 
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where 

b  *  -  [r^  +  vKS-j+r^l  [&jft2  +  ^3^ 

-  y r2  [ 2«1a>2  -  O^2  +  $2\  + 

and 

c  *  -  [r3  +  y(A^+r2)]  [^2^3  "  ^2^ 

-  ur2[2^1a)3  +  Q^2  +  0^  -  u2]  , 

Equation  for  Angular  Momentum 

H  -  Hx  ix  +  H2  Jx  +  H3  gx  ,  (D-31) 

where 

H1  *  (Ai+Bl)fil  *  (D-32a) 

H2  “  '03[B2  +  °  r2  a1+r2)“Bl]fJl  +  [B2  +  °  r2(Vr2)la)2 

+  [A2  +  B2  +  a(i1+r2)2]fi2  (D-32b) 

and 

H3  ”  ©2CB2  +  a  r2(£1+r2)-B1]fi1 

+  [B2  +  a  r2(2.1+r2)]u>3  +  [A2  +  B2  +  a^l+r2^S  .  (D-32c) 

Equations  for  the  Angles  9^  and  y 

*  Q  +  02  cos*  -  03  sin*  (D-33a) 

and 

fN  -  4*  +  02  sin*  +  03  cos*  .  (D-33b) 
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